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FOREWORD 


This  user's  manual  was  prepared  by  the  General  American  Research  Division 
of  the  General  American  Transportation  Corporation  for  the  Stanford  Research 
Institute  under  Subcontract  No.  11599 (63OQA-O9O) •  It  explains  in  detail  the 
procedure  for  using  the  ‘Mr. outer  code  to  calculate  the  response  of  a  shelter 
to  time -varying  thermal  loads  and  time-varying  ventilacion  inlet  psychrcmetric 
conditions.  The  code  is  based  upon  research  performed  between  June  1963  and 
May  1968  on  the  Analysis  of  Shelter  Ventilation  Requirements  which  falls  in  the 
program  area  of  OCD  Work  Unit  1215A . 

The  project  was  monitored  by  Mr.  F.  Allen  (then  rt  OCD)  and  Mr.  D.  Bettge 
of  OCD,  and  Messrs.  C.  Grubband  J.  Halsey  of  SRI. 
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and  A.  Kapil  of  GARD,  Dr.  J.  Buchanan  of  OCD,  Mr.  P.  Achenbach  and  Dr.  T.  Kusuda 
of  the  National  Bureau  of  Standards,  Mr.  W.  Spiegel,  Consulting  Engineer,  and 
Messrs.  F.  Hughes-«.aley  and  T.  Hori  of  SRI  for  their  comments  and  criticism. 

The  valuable  cooperation  of  the  personnel  at  the  OCD  computer  facility  in 
Olncy,  Maryland  and  the  Control  Data  Corporation  computer  facility  in  Chicago, 
Illinois  is  acknowledged. 

It  should  be  noted  that  the  version  of  the  5KEP  code  presented  in  this 
report  (Mod.  3)  has  been  superseded  by  a  new  modification  (Mod.  6).  This  new 
code  employs  a  slightly  different  sequence  of  calculation  to  dacruase  calcu¬ 
lation  time. 
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ABSTRACT 


The  Shelter  Environmental  Prediction  (SlfEP)  Computer  Code,  Modification  Z, 
calculates  the  response  of  a  shelter  to  time-varying  thermal  loads  and  time- 
varying  ventilation  inlet  psychrometric  conditions.  The  code  accommodates  solar 
radiation,  boundary  surface  heat  transfer,  moisture  condensation  on  the  boundaries, 
equipment  and  lighting  loads,  and  air  conditioning.  The  manual  explains  the 
various  features  of  the  code,  presents  the  input  and  output  formats,  gives 
representative  input  data  values,  and  includes  an  application  of  the  code  to 
a  typical  shelter  structure. 

The  required  inputs  to  the  code  and  the  outputs  obtainable  from  it  are 
given  belov; 


Inputs 


1. 

Ventilation  rate 

1. 

hourly  temperatures 

2. 

Inlet  psychrometric  con¬ 
ditions 

2. 

Hourly  air  conditioning  load 

3- 

Shelter  geometric  charac¬ 
teristics 

3. 

Hourly  boundary  heat  transfer 

4. 

Shelter  thermal  charac¬ 
teristics 

4. 

Hourly  metabolic  loads 

5. 

Shelter  latitude,  longitude 
and  altitude 

5. 

Average  temperatures  (dry-bulb,  vet. 
bulb,  effective)  for  the  entire  , 
occupant  or  specified  hours  of 
occupancy 

6. 

Period  of  occupancy 

6. 

Maximum  and  minimum  temperatures 
(dry-bulb,  wet-bulb,  effective) 
during  the  occupancy. 

The  '•ode  it  written  in  FORTRAN  Vf  specifically  for  the  CSC  36 00  coeputer 


system  and  consists  of  approximately  ?60  cards.  It  requires  a  memory  bank  of 
et  least  65*  words  (each  word  of  48  bits)  for  the  code  and  the  accompanying 
matrices.  For  the  CDC  3600  erd  for  a  problem  comparable  to  the  sample  problem 
presented  in  the  text,  the  compilation  time  is  approximately  1  min.,  4$  sec. 
end  the  execution  time  approximately  1  sec.  for  each  hour  of  shelter  occupancy. 
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SECTION  1 
INTRODUCTION 

The  Shelter  Environmental  Prediction  (SKEP)  code  will  calculate  the 
environmental  response  of  a  shelter.,  or  any  large  structure,  to  time-varying 
thermal  leads  and  ventilation  inlet  psychrometric  conditions.  Shelter  con¬ 
ditions  computed  by  finite  difference  methods  define  the  ui story  of  the  shelter 
environment  dui’ing  occupancy. 

The  SHEP  code  accommodates  solar  radiation,  boundary  surface  heat  transfer, 
moisture  evaporation  or  condensation  on  the  boundaries,  equipment  and  lighting 
loads,  and  air  conditioning  .{derivations  of  the  corresponding  equations  are 
given  in  the  Appendices).  Little  technical  knowledge  of  mass  or  energy  trans¬ 
fer  is  required  bv  a  user  of  the  SHEP  code.  Information  that  must  be  supplied 
includes  shelter  physical,  thermal  and  geometrical  characteristics,  occupancy 
levels  and  durations,  and  hourly  inlet  conditions.  The  amount  of  required  data 
is  admittedly  large,  hut  is  necessary  so  that  the  SHEP  code  will  remain  appli¬ 
cable  to  a  vide  variety  of  shelters. 

The  shelter  a 'del  used  in  the  code  is  based  upon  the  assumptions  that: 

1.  Air  within  thn  shelter  end  ventilation  air  is  completely  and 
instantaneously  mixed,  to  that  one  psychrometric  erudition  will 
completely  specify  the  shelter  environment} 

1.  Moisture  condensation  and  evaporation  is  filmvit*  with  the  con¬ 
vective  heat  transfer  coefficients  constant  over  etch  boundary 
surface,  but  variable  from  boundary  to  boundary. 

3.  Radiative  energy  transfer  between  the  shelter  boundaries  is 
neglected. 

BCWCRAl  AM e*IIC  AM  KCttAMCM  DtVHHOW 
1-1 


4. 


Shelter  exhaust  air  is  at  the  psychrometric  condition  of  the 
shelter  atmosphere;  and 

5.  Thermal  and  physical  properties  of  the  air  and  of  the  structural 
nate rials  are  not  temperature -dependent. 

The  SHEP  code  permits  the  analysis  of  an  actual  shelter  vith  up  to  20  different 
boundaries,  enclosing  a  volume  which  has  sensible  and  latent  energ*-  exchange 
with  the  environment  external  to  the  shelter. 

This  user's  manual  explains  the  various  features  of  the  SHEP  code,  presents 
the  input  and  output  formats,  gives  representative  input  data  values,  and  in¬ 
cludes  an  application  of  the  code  to  a  typical  shelter  structure. 

The  code  is  written  in  FORTRAN  IV  specifically  for  the  CDC  36OO  computer 
system  and  consists  of  opproxiaately  780  cards.  It  requires  a  memory  bank  of 
at  least  65*  words  (each  word  of  48  bits)  for  the  code  and  the  accompaiying 
matrices. 

Pv  the  CDC  36CO  and  for  a  problem  comparable  to  tbe  sanq>le  problem 
presented  in  the  text,  the  compilation  time  is  approximately  1  sin.,  45  sec. 
and  the  execution  time  approximately  1  sec.  for  each  hour  of  shelter  occupancy. 
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SECTION  2 

FEATURES  OF  THE  PROGRAM 

The  shelter  model  used  by  the  SHEP  code  is  an  accurate  description  of  the 
actual  shelter  under  analysis  because  of  the  methods  which  have  been  formulated 
in  the  code  to  handle  the  thermal  processes  which  affect  the  shelter  environment. 
These  processes  include: 

1.  solar  radiation, 

2.  air  conditioning, 

3.  moisture  evaporation  and  condensation, 

4.  time-varying  inlet  conditions,  and 
5*  boundary  heat  transfer. 

It  has  been  found  that  during  extreme  hot  weather,  many  actual  shelters  respond 
as  if  they  had  adiabatic  boundaries.  Thus,  the  code  has  the  ability  of  analyzing 
a  shelter  with  adiabatic  or  non-adiabatic  shelter  boundaries. 

The  complexity  of  the  code  has  necessitated  the  printout  of  diagnostic 
error  messages.  These  are  included  as  a  safeguard  against  using  illegal 
variable  values,  or  exceeding  the  inherent  limits  of  thvt  equations  comprising 
the  code . 

2.1  Shelter  Boundaries 

An  actual  shelter  it  bounded  by  walls,  floor  end  ceiling,  composed  of 
several  different  kinds  of  material*  (e.g.*  cement,  gl»*.«,  metal,  etc.)  and 
these  can  be  at  any  orientation  (tJope  and  azimuth).  In  addition,  the  exterior 
boundary  surface  may  be  expceed  to  various  media  (e.g.,  elr,  soil,  etc.), 
therefore,  a  "boundary”  is  dafined  as  a  plane  or  approximately  plane  curved 
surfre.e*  of  the  shelter  enclosure  consisting  of  e  homogeneous  inner  layer  a 

•See  Appendix  E. 
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successive  homogeneous  layers  exposed  to  one  external  medium.  If  windows  are 
included  in  one  wall  of  the  shelter  and  this  wall  is  partially  underground, 
then  this  wall  consists  of  three  boundaries  (see  Figure  2-1).  A  shelter  volume 
may  be  defined  by  up  to  20  boundaries.  A  curvilinear  boundary  should  be  con¬ 
sidered  plane,  with  its  orientation  that  of  the  plane  tangent  to  the  boundary 
at  its  midpoint  (see  Appendix  s). 

An  implicit  finite-difference  approach  based  upon  a  one -dimensional  node 
array  is  used  to  compute  temperature  distribution  across  a  boundary,  with  each 
boundary  spanned  by  up  to  ho  nodal  points.  A  boundary  may  consist  of  up  to 
5  layers,  each  having  distinct  thermal  properties  and  nodal-point  spacings. 

An  air  space  within  a  boundary  and-a  soil  media  adjacent  to  a  belowgrade 
boundary  are  considered  &3  layers  of  the  respective  boundaries.  Shelter  tests 
have  shown  that  10  feet  from  the  exterior  boundary  surface,  the  soil  temperature 
is  invariant.  Thus  the  soil  can  be  considered  as  a  layer  10  feet  thick  whose 
exterior  surface  is  at  the  undisturbed  soil  temperature. 

Nodtl  point  spacing  in  a  typical  boundary  consisting  of  k  layers  is  shown 
in  Figure  2-2.  Each  layer  is  divided  into  slab  A thick  end  a  nodal  point  is 
placed  at  the  center  of  each  slab.  The  Innermost  end  outermost  points  are 
Ax^/2  from  the  next  inner  nodal  point  and  are  placed  at  the  surface  of  the 
layer  (see  Appendix  C).  It  should  be  noted  that  thare  is  e  nodal  point  at 
each  surface  and  at  <t«ch  interface,  and  that  the  nodal  point  numbering  is 
continuous . 

The  number  of  nodal  points  in  each  layer  is 

*  .  ;j*»r -thickness  +  A 

»  A  X 

k 

which  includes  the  node  on  the  inner  surface  of  the  layer.  The  number  of 
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nodal  points  in  a  wall  excluding  the  outermost  point  is  HMAX.  and  must  be 

V 

supplied  by  the  programmer.  It  is  given  by 


UMAX, 


Ls 


layers 


and  the  total  number  of  nodal  points  is 


H  -  UMAX.  +  1 

J  V 

Spacing  between  nodal  points  should  vary  inversely  with  conductivity. 
Suggested  spacing  in  soil  is  one-haJf  to  one  foot;  in  structural  materials, 
one  inch. 

The  nodal  point  approach  to  the  temperature  profile  across  the  boundary 
permits  an  implicit  representation  of  the  temperature  through  the  boundary 
via  the  simultaneous  linear  equations: 

•llTl  +  *12T2 
*21T1  4  *22T2  t  *23T3 

*32T2  4  *33T3  +  *3UTU 


*n-l,n-2Tn*2  4  Vl,a-lTn-l  4 
*n,n-l*n-l 

obtained  by  applying  an  energy  balance  ot  each  of 
the  coefficients  the  temperatures  at  each  node 
storage  terms,  which  are  ties -dependent.  In  matrix  form,  the  equations 
become: 


*n-l,n*n  ^o«l 

*mt„  ■  \ 

n  notes.  The  a*»  ere 
sod  the  b*s  are  the  heat 
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These  equations  may  be  solved  for  the  nodal  temperatures: 

M  -  W'1  [»]  • 

and  the  heat  transfer  through  a  non-adiabatic  boundary  due  to  conduction  moy 
be  computed  based  on  the  interior  boundary  surface  temperatures.  In  addition, 
for  non-adiabatic  boundaries,  heat  loss  due  to  moisture  condensation  or  evapor¬ 
ation  on  the  interior  boundary  surfaces  is  computed.  The  derivations  of  the 
equations  for  determining  the  £aJ  and  jfij coefficients,  and  those  for  heat  trans¬ 
fer  due  to  condensation,  are  given  in  Appendix  C.  This  implicit  procedure  does 
not  have  the  stability  considerations  associated  with  the  explicit  method. 

Thus,  longer  time  intervals  can  be  utilised  at  each  finite  difference  inc remeat. 
If  the  boundaries  are  considered  adiabatic,  heat  transfer  both  through  the 
boundary  and  due  to  condensation  or  evaporation  on  the  surface  are  considered 
zero. 

2.2  Solar  Radiation 

..  mmmm mmmm  «m  m  —  ■■  i 

If  a  surface  is  exposed  to  ambient  air,  It  will  normally  be  exposed  to 
direct,  diffuse  and/or  reflected  sola:  radiation.  If  the  surface  ia  a  shelter 
boundary,  the  incident  radiation  will  serve  to  raise  the  shelter  temperature 
in  two  ways.  First,  '.uuident  radiation  will  be  psrtlsUy  absorbed  by  the 
boundary,  thus  raising  its  temperature  and,  In  time,  the  temperature  of  the 
shelter.  Second,  radiation  transmitted  through  the  boundary  will  be  incident 
on  interior  surfaces,  similarly  raising  thsir  temperatures  and  that  of  the 
ahelter  atmosphere. 

SHE?  makes  several  simplifying  assumptions  in  treating  the  effects  of 
solar  radiation.  One  of  these  ia  that  all.  reflected  radiation,  both  that 
which  fall*  onto  the  ahelter  surfaces  from  surrounding  objects  end  that  which 
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is  reflected  from  the  shelter  surfaces,  is  neglected.  Also,  radiation  trans¬ 
mitted  through  a  boundary  is  considered  an  instantaneous  thermal  load  inside 
the  shelter,  similar  to  the  lighting  and  equipment  loads. 

All  calculations  involving  solar  radiation  are  done  using  solar  time, 
which  differs  slightly  from  Civil  Time.  The  relation  may  be  expressed  as 

8olar  Time  *  Civil  Time  +  Equation  of  Time  +Ajj0ng 

Civil  Time  is  that  time  assigned  to  a  zone  covering  approximately  15°  of 
longitude;  it  is  the  same  as  Standard  Time.  Since  local  time  actually  does 
not  remain  constant  within  each  15*  tone,  but  varies  by  four  minutes  for  each 
degree  within  the  tone,  a  correction  factor, Aj ^  t  takes  this  into  account. 

The  Equation  of  Time  is  the  difference  between  mean  solar  time  and  apparent 
solar  time  and  is  given  ss  the  (l^-Epheoeris  Transit).  See  Reference  1. 

Solar  time  is  computed  in  tha  program  knowing  the  time  of  the  start  of 
occupancy,  the  longitude  of  the  shelter,  and  tha  first  day  of  occupancy. 

Variation  of  the  Equation  of  Tima  is  shown  in  Figures  2-3  »nd  the  values 
used  in  the  program  are  listed  in  Table  1.  The  daily  veitiee  of  the  Equation 
of  Time  chenge  alightly  from  year  to  year,  but  the  actual  values  remain  within 
0-3  min.  of  those  listed  (see  Reference  2). 

Solar  rsdistlon  intensity  is  computed  ones  each  time  increment  tor  tach 
exposed  wall  (See  Appendix  B.)  Although  radiation  will  penetrate  into  * 
surface  a  certain  distinct  before  being  absorbed,  the  SHIP  code  assumes  this 
absorption  phenomenon  is  s  surface  effect.  Also,  even  though  the  transmittance 
of  the  boundary  eater) al  may  vary  *ltl»  incident  waveler,gth,  only  a  total  trans¬ 
mittance  value  is  used  in  the  code. 
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2.3  Air  Conditioning 

In  »  typical  ftir  conditioning  unit,  warm  ftir  passer  over  cooling  coils, 
thereby  lowering  its  dry-bulb  temperature  through  convective  eat  exchange 
and  lowering  its  moisture  content  through  condensation.  TUo  temperature  of 
the  coils  is  not  constant,  but  varies  ou  a  function  of  load.  Typically,  all 
the  air  passing  through  the  air  conditioner  does  not  experience  these  heat  and 
mass  transfer  processes  due  to  flow  turbulence. 

The  method  the  SHEP  code  utilises  to  calculate  the  air  conditioning  load 
assumes  that  the  temperature  of  the  coils  and  the  percentage  of  air  unaffected 
by  passing  through  the  unit  do  not  vary  with  the  cooling  load.  An  "effective 
coil  surface  temperature"  ft  assumed  as  the  design  temperature  of  the  coils; 
and  a  "coil  bypass  factor"  is  Use  design  value  of  the  percentage  of  air  un¬ 
affected  by  the  cooling  unite.  If  theee  two  parameters  are  both  assumed  con¬ 
stant,  the  air  conditioning  performance  of  an  actual  unit  can  be  estimated 
reasonably  well-  See  Appendix  F,  which  is  based  upon  Reference  3* 

A  typical  air  conditioning  system  for  s  shelter  may  be  represented  by 
the  schematic  diagram  of  Figure  2-**;  the  various  psychrometric  states  of  tho 
air  in  the  system  are  indicated  in  the  accompanying  psychrometric  chart. 

Rote  that  "return  air  bypass"  refers  to  air  not  pasting  through  air  condition¬ 
ing  ur.lt;  "coil  bypass  factor”  la  a  characteristic  of  the  unit  itself.  Outside 
air  at  state  point  1  enters  the  system  end  is  adiabstlcally  mixed  with  return 
air  at  the  shelter  condition,  3,  to  create  state  point  2.  A  portion  of  this 
air  peases  through  the  air  conditioning  unit,  exiting  with  e  lower  dry-bulb 
temper* tuft  and  «  low  humidity  ratio,  state  point  3,  end  is  then  adiabetlcally 
mixed  with  the  unconditioned  air  to  give  state  point  b.  Air  at  chit  conditio* 
is  the  inlet  air  to  the  shelter  which  is  then  used  to  compute  the  new  shelter 
condition,  S'. 
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The  ventilation  rates  of  outside  ait**entering  the  system,  of  air  passing 
through  the  air  conditioning  unit,  and  of  air  entering  the  shelter  must  all 
be  specified  by  the  code  user.  Therefore,  the  mixing  processes  can  be  altered, 
or  even  eliminated,  by  varying  the  various  ventilation  rates. 

Air  conditioning  removes  both  sensible  and  latent  energy.  This  energy 
removal,  AQ,  is  given  by 

AQ  *  c  mAT  +  h-  mAW 
P  fg 

where 

Cp  *  specific  heat  of  air«  0.24  Btu/lb  °P 

m  *  mass  flow  rate  of  dry  air  through  the  air  conditioning 

unit,  lb /hr 

am  * 

AT  ■  temperature  difference  across  the  unit,  *F 

h^  *  heat  of  condensation  of  water  vapor,  Btu/lb  vater 

aV  -  humidity  ratio  difference  across  the  unit,  lb  water /lb  dry  air 

2.4  Moisture  Condensation 

Wherever  the  vapor  pressure  at  the  interior  surface  of  a  boundary  is 
lover  than  the  vapor  pressure  of  the  shelter  atmosphere,  moisture  will  migrate 
to  and  condense  on  the  surface,  thereby  removing  energy  from  the  shelter  envir¬ 
onment.  If,  at  another  time,  the  vapor  pressure  at  the  boundary  surface  is 
higher  then  the  vapor  pressure  of  the  shelter  air,  the  proceaa  will  be  reversed 
and  condensed  water  on  the  surface  will  re-tveporete  and  migrate  into  the 
shelter  atmosphere. 

If  sufficient  moisture  condenses,  it  will  ran  down  the  wells  end  accumu¬ 
late  on  th|  floor  of  the  shelter.  Since  the  SHEP  code  assume'*  that  condensation 
occur*  in  an  even  film  on  an  interior  surface,  the  effect  of  grav'ty  is 
neglect* i. 
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In  the  code,  the  vapor  pressure  at  the  boundary  surface  is  assumed  to  be 
the  saturation  pressure  of  water  vapor  at  the  temperature  of  the  boundary 
surface.  Values  of  saturation  pressures  as  a  function  of  dry -bulb  temperature 
(see  Table  II )  are  a  part  of  the  standard  input  data.  Only  the  range  of  tem¬ 
peratures  from  50*F  to  120*F  are  considered,  since  metabolic  relations  used  in 
the  SHEP  code  are  invalid  outside  of  this  range. 

2.5  Metabolic  Loads 

The  metabolic  expressions  of  Hough ten?  are  adopted  in  this  study  since  they 
reflect  effects  of  relative  humidity  in  the  shelter.  They  are  given  by 


Wbl.  “  -0.06875(Tdb)2  .  1.625(Tdb)  .  523.0 

I-1.h62  (et)  ♦  Bl^.O 
-1.508  (ET)2  4  259.7  (ET)  -  10795.2 
0.0  (assumed) 


for  T..  >  50*F 

CD  *- 

for  50*F  <  ET  <  87*F 
for  8?*F  <  ET  <  102 *F 
for  ET  >  102  *F 


vith  latent  *  Qtotal  '  Sensible*  th#t  tke  bo^  is  ir  thcrfr#1  e*uilib' 

rium  with  the  shelter  environment,  and  hence,  neglecting  the  heat  storage  term. 

2.6  Auxiliary  I.oads 

Lighting  and  equipment  loads  tre  introduced  ss  inputs  and  may  be  tine- 
varying  (on  *r,  hourly  basis). 

2.7  Varying  Inlet  Conditions 

The  psych rometric  condition  of  the  aheltcr  atmosphere  and  of  the  entering 
ventilation  air  are  considered  by  the  code  as  a  function  of  dry-buib  temperature 
and  relative  humidity. 

Certain  factors  which  effect  the  shelter  atmosphere  normally  very  with 


time;  art  lent  temperature  and  relative  tmidity  ara  chief  among  these.  Other 
far?  ?r:  vM'-h  :.ay  charge  during  occupancy  arc: 


’Superscripts  refer  to  references,  p.  7-1. 
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TABUS  II 


SATURATION  PRESSURE,  Pg,  AS  A  FUNCTION 
OF  DRY -BULB  TEMPERATURE,  Tdb 


Tdb* 

P  *K 

rs 

T 

ab 

PS 

Tdb 

P8 

T 

Adb 

ps 

50 

25-648 

64 

42.494 

78 

68.357 

92 

107.021 

106 

163.368 

51 

26.620 

65 

44.006 

79 

70.646 

93 

110.390 

107 

168.235 

52 

27.6T2 

66 

45.562 

80 

72.994 

94 

113.846 

108 

173-218 

53 

28.656 

67 

47.174 

81 

75.413 

95 

117.403 

109. 

178.330 

54 

29-724 

68 

48.816 

82 

77-904 

96 

121.061 

110 

183.571 

55 

30.830 

69 

50.530 

83 

Co.  467 

97 

124.819 

111 

188.942 

56 

31.968 

70 

52.286* 

84 

"83.102 

98  * 

128.664 

112 

194.458 

57 

33.149 

71 

54.086 

85 

85.810 

99 

132-624 

113 

200.102 

58 

34.358 

72 

55-958 

86 

88.589 

100 

136.685 

114 

205 .891 

59 

35.611 

73 

57.873 

87 

91.454 

101 

140.846 

115 

211.809 

60 

36. 90; 

74 

59-846 

88 

94.406 

102 

145.123 

116 

217.672 

61 

38.232 

75 

61.891 

89 

97.430 

103 

149.501 

117 

224.107 

C2 

39.6J4 

76 

63.979 

90 

100.541 

104 

154.008 

118 

230.486 

63 

41.040 

77 

66.139 

91 

103.738 

105 

158.630 

119 

237. 010 

120 

243 .706 

*»  ,  5 

p§:  ib/rr 
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1.  the  number  of  occupants, 

2.  the  ventilation  rate,  and' 

3.  the  lighting  and  equipment  loads. 

In  the  code,  these  factors  may  have  different  values  for  each  hour  of  occupancy. 

2.8  Error  Messages 

Since  mar.y  of  the  relations  used  in  SHEP  are  valid  only  over  certain  ranges, 
limits  have  been  placed  on  the  values  of  certain  input  data  and  of  certain  com¬ 
puted  variables.  If  these  limits  are  exceeded,  diagnostic  error  messages  are 
printed,  stating  what  variable  is  in  error  and  vhat  value  it  has  at  the  time 
it  exceeded  its  limit.  If  an  input  variable  has  an  illegal  value,  this  will  be 
indicated  and  the  run  will  terminate ;“if  a  computed  variable  la  invalid*,  this 
fact  will  be  printed  and  the  run  will  continue.  A  third  type  of  diagnostic 
message  is  generated  through  machine  error;  lor  instance,  if  an  index  has  over¬ 
run  its  bound. 

A  list  of  the  error  messages  is  indicated  in  Table  III.  These  error 
messages  arc  incorporated  into  tht  SHEP  code  and  are  in  addition  to  the  com¬ 
pilation  and  execution  diagnostic  routines  of  the  machine. 

2 . 9  Petcrr  lrs "ion  of  Shelter  Dry-Bulb  Tetnper.v,uro 

All  titer;  al  loads  and  enthalpies  arc  computed  at  the  begining  of  a  time 
increment  based  upon  an  estimate  of  what  the  shelter  dry-bulb  temperature  will 
be.  Then  the  net  amount  of  energy  added  to  the  shelter  system  during  that 
increment  is  calculated.  With  this  value,  the  shelter  dry -bulb  temperature 
can  be  computed  by  equation  20  of  Appendix  A.  Since  the  shelter  dry-bulb 
tespore’urc  of  the  inerervent  is  initially  estimated  to  be  that  of  the  last 
tner^ent,  a-,  'tereticn  procedure  is  used  to  permit  each  calculated  temperature 


CfiVfiiM.  AwcnicAh  nr.ir.Ancn  mvortov 


2-15 


TABLE  III 
ERROR  MESSAGES 

1.  Input  Variables  Illegal 

A.  "Inclination  Angle  XI  ( _ )  ic  _ .  This  value  muat  be  90°  or  leBS." 

B.  "The  latitude  of  the  shelter  is  greater  than  the  program  can  handle 
for  this  day  of  the  year." 

2.  Calculated  Variables  Illegal 

A.  "R  during  hour _ is  _ ,  which  is  less  than  1.0.  Therefore,  the 

latent  metabolic  energy  has  been  reduced." 

B.  "Shelter  effective  temperature  during  hour  _  is  _  which  exceeds 

the  limits  on  the  metabolic  relations.1' 

C.  "Shelter  dry-bulb  temperature  during  hour  _  is  _  which  exceeds 

the  limits  on  the  metabolic  relations." 

D.  "The  relative  humidity  in  the  shelter  has  exceeded  saturation." 

3*  Procedural  Errors 

A.  "The  limits  on  layaax  or  nmax  of  wall  _  have  been  exceeded." 

"The  number  of  shelter  boundaries,  NW,  has  been  exceeded.” 

C.  "The  number  of  iterations  for  an  increment  has  exceeded  the  value 

cf  IUATA(3)." 

P.  "The  number  of  increments  in  an  hour  exceeds  the  value  of  IDATA(5)." 

B.  "The  hour  numbers  are  not  indexing  properly." 

K.  "The  hour  number  has  exceeded  the  length  of  occupancy." 
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to  be  used  as  a  new  estimate  to  compute  a  more  accurate  shelter  dry-bulb 
temperature . 

The  number  of  Iterations  for  each  increment  is  a  function  of  the  accuracy 
of  the  initial  estimate.  If  the  dry -bulb  temperatures  computed  on  two  succes¬ 
sive  iterations  are  less  than  l’F  apart,  the  last  temperature  computed  is 
assumed  to  be  the  value  for  the  increment  and  the  iteration  procedure  is  cur¬ 
tailed.  Otherwise,  the  procedure  is  continued  until  the  number  of  iterations 
reaches  the  value  of  IDATA(3),  (see  Section  3)* 

2.10  Determination  of  Shelter  Effective  Temperature 

After  the  shelter  dry-bulb  temperature  for  an  increment  has  beer,  computed, 
wet-bulb  temperature  is  determined  via  the  Carrier  equation  (equation  4l  of 
Appendix  A).  Using  a  linear  approximation  to  the  nomogram  of  Reference  4,  the 
effective  temperature  is  calculated  from: 


ET 


107.5  Tdb  -  1.5.2  Tub  j  1.5'F  <  Tab  <  110‘F 

for 


Tdb  -  Twb  +  62‘3 


30°F  <  Twb  <  110*F 


where  T *  ary-buib  temperature,  *F 
00 

*  vet-bu.lb  temperature,  *F 

The  psychrometric  condition  of  the  shelter  atmosphere  at  the  beginning  of 
the  next  increment  is  then  set  equal  to  that  at  the  end  of  the  lest  increment 
and  calculations  for  the  next  increment  are  begun.  When  one  hour's  calcu¬ 
lations  have  been  completed,  the  values  for  the  last  increment  of  the  hour  arc 
printed  out  as  the  condi tions  at  the  end  of  the  hour. 
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SECTION  3 
INPUT  DATA 


Input  data  consist  of  physical  constants,  hourly  occupancy  constants, 
geometvicai  and  thermal  characteristics  of  the  shelter,  and  initial  conditions. 
The  input  variables  and  their  definitions  are  given  below;  where  units  are  not 
indicated,  the  value  is  dimensionless. 


ALAT 

ALONG 


-Latitude  of  the  shelter,  degrees.  The  limits  on  th«  range  of 
ALAT  are  a  function  of  the  time  of  the  year,  see  Figure  3-1. 


-Longitude  of  the  shelter,  degrees. 
0°  to  l8o#  W  ere  allowed. 


5rly  values  in  the  rswge 


ALT 

AREA(j) 

ASOIL(j) 


-Altitude  of  the  shelter,  feet. 

2 

-Area  of  boundary  -j,  ft  .  - 

•Control  constant  determining  the  media  exterior  to  the  outer 
surface  of  boundary 


ASOIL(j) 


(-1  arableni 
0  inte  r' 
1  oi* 


-1  ambient  air 
or  air 


-Arimuth  angle  of  boundary  * .  AZ(j)  ia  measured  in  a  clockwise 
direction  from  north  to  the  out ward -pointing  normal  to  the 
boundary.  Values  between  0*  and  360*  are  permitted.  For  a 
boundary  facing  east,  A7.(j)  *  30*;  if  it  facet  vest,  A Z($)  * 
270*. 


•Ceding  coll  bypass  factor  of  the  air  conditioning  unit.  It  is 
defined  as 


8F 


where 


.  iLLi 

‘•d  - 

•  air  temperature  leering  the  coil,  *F 

t  *  coil  effective  surface  temperature,  *F 
*  (TCOIL) 

t  *  air  temperature  entering  the  coil,  *F 


The  suggested  vdue  ia  0.1V- 


nrMcnAu  AfwcAfj  nnsrAicH  orvivatw 
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Figure  3-1  Allowable  Values  of  the  Shelter  Latitude 


CfKTN 


COK(k,j) 

CP(k,j) 

D 

DP(k,j) 

DT 

DX(k,j) 

EQTlM(i) 

PI 

P2 

GG 

HK 

Hl(j) 

HO(j) 

HRS 

IDATA(l) 

It*TA{2) 

I»TA(3) 

1MTA(*0 


-Total  ventilation  rate  through  air  conditioner,  cfm. 

-Thermal  conductivity  of  material  of  layer  k  of  boundary  j. 
Btu/hr-ft-*F. 

-Specific  heat  of  material  of  layer  k  of  boundary  j,  Btu/lb-*P. 

-Day  number  of  the  start  of  occupancy.  Values  in  the  range  of  1-36$  are 
allowed,  e.g.,  for  January  1,  D  ■  1. 

-Thickness  of  layer  k  of  boundary  j,  feet. 

-Time  increment  of  shelter  history  used  in  finite-difference 
equations,  hrs.  0.00  <  DT  <  0.500  hours,  DT  *  0.16?  is 
suggested.  ~ 

-Nodal  point  spacing  of  layer  k  of  boundary  J,  feet. 

-Equation  of  Time  for  day  i,  hours.  See  Table  I. 

-Total  'entilation  rate  entering  shelter,  cfm. 

-Total  ventilation  rate  of  outside  air  entering  shelter  system,  cfm. 

-Equivalent  into, ior  lengti  of  shelter,  feet.  Equivalent  dimensions 
are  those  of  a  parallelepiped  approximating  the  shelter  under 
analysis. 

-Equivalent  interior  height  of  shelter,  feet. 

^Interior  film  heat  transfer  coefficient  of  boundary  j, 

Btu/hr-ft  -*F. 

•Exterior  film  heat  transfer  coefficient  of  boundary  j, 

Btu/hr-ft*-*F. 


-Total  number  of  continuous  hours  of  shelter  occupancy. 


•Sets  initial  boundary  temperatures 


2I*TA(1) 


1  -  all  boundary  nodes  have  different 
temperatures 

0  -  all  boundary  nodes  have  same 
temperature 


-Reference  constant;  IDAT*(2)  *  Os  suggested. 


-Kvcier  of  iteiationa  for  computing  incremental  shelter  dry-bulb 
temperature.  I r*TA ( 3 )  -  is  suggested. 

•Srier  of  shelters  to  be  snalyred  in  the  cccpui-ir  run. 


SENriVI.  XMEmCAN  IKlI-iU  N 
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4  • 

jATA(5)  -1  hour/DT,  rounded  to  the  nearest  integer. 


IDATA(b) 


-Sets  type  of  boundary . 


IDATAvo) 


for  all  non-adiabatic  boundaries 
for  all  adiabatic  boundaries 


IDATA(?) 


-Determines  presence  of  air  conditioning. 


IDATA(7) 


0  -  no  air  conditioning 
1  -  air  conditioning  is  present 


XDATA(8)  -Reference  constant,  It.aTA(8)  *  0  is  suggested. 


IDATA(9)  -Reference  constant,  BATA(<y)  *  0  is  suggested. 


IDATA(IO)  -Sets  atmospheric  condition  for  solar  radiation. 


IEATA(IC)  * 


6 


-  for  clear  atmosphere 

-  for  industrial  atmosphere 


IOaXA(ll)  -Establishes  format  of  printout. 


IDATA(Ii)  * 


1  -  Standard  output  format 
'  X)  -  Only  ftaximibn,  minimum  and  average 
„  temperatures  ore  printed 


IDATA(12)  -Establishes  variab'Uty  of  occupancy  values. 

IDATA(12)  ■  CilGS,  Q1ITE,  QEQUIP  change  for  each  hou: 


These  variables  are  held  constant  ror  entire 
occupancy. 


IUATA(13)  -Establishes  variability  of  Inlet  values. 


IMTA(13) 


Inlet  condition#  held  constant  for  entire 
occupancy 

These  variables  change  for  e.ary  hour  of 
occupancy. 


IET(m)  -Hour  lie.it*  of  the  Interv&l*  for  which  average  temperatures  ate 

desired,  e.g..  for  average  dry-bulb  temperatures  from  hour  number  3 
to  hour  number  9  of  occupancy,  then  llft(l)  *  3,  IHT{2)  *  9.  IKT(n), 
whore  »  is  an  odd  integer,  1  <  o  <  5?  determine*  the  atari  of  an 
interval;  and  1KT(«*1)  deUn«Tne*~the  end  of  the  interval. 


LA7KAX(j)  -Kunber  of  layers  in  boundary  j.  .toil  back  of  a  boundary  or  an 

air  space  within  a  boundary  is  considered  a  layer  of  the  boundary. 
Value  must  be  5  or  less. 


U«T(1)  -local  Standard  Time  at  beginning  of  occupancy,;  defined  in  terns 
of  ailiUry  time  (0-23  hrs). 


OCNPniM.  AW  mcAM  RES  'M.CK  bivt&io^ 
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NMAX(j)  -  Number  of  nodal  points  in  boundary  j.  If  a  boundry  is  soil- 
backed,  the  nodal  points  extend  into  and  through  the  soil,  up 
to  a  suggested  dista -oe  of  30  feet.  Xntegrel  values  up  to  39 
are  permitted.  The  outermost  nodal  point  is  not  included  in 
this  value. 

NOINT  -  Number  of  time  intervals  for  which  average  temperatures  are 
desired,  NOINT  <  30. 

NW  -  Number  of  boundaries  in  the  shelter,  KW  <  SO. 

P  -  Number  of  occupants. 

A 

PSAT(n)  -  Saturation  pressure  of  water  vapor  at  n  *F,  lb/ ft  . 

See  Table  II. 

QBQUIP  -  Sensible  heat  load  due  to  mechanical  and/or  electrical  equipment 
(excluding  lighting),  Btu/hr. 

QLITE  -  Sensible  heat  load  due  to  lighting,  Btu/hr. 

RCW(k,j)  -  Density  of  mate_rial_of  layer  k  of  boundary  J,  lb/ft  . 

Kll  -  Relative  humidity  of  inlet  air  at  beginning  of  occupancy. 

R12  -  Relative  humidity  of  inlet  air  at  end  of  each  hour  of  occupancy. 

R2  -  Relative  humidity  in  shelter  at  beginning  of  occupancy. 

TA1  -  Dry-bulb  temperature  of  air  exterior  to  shelter  structure  at 

beginning  of  occupancy,  *F. 

TA2  -  Dry-bulb  temperature  of  air  exterior  to  shelter  struc*/.re  »t 

end  of  each  hour  of  occupancy,  *F. 

TCOIL  -  Sffective  coil  temperature  of  air  condi tioning  unit,  ’F. 

Suggested  value  it  5>*F  vith  VCOIL  *  0.00921  'b  water/lb  dry  air. 

TD11  -  Dry-bulb  temperature  of  inlet  air  at  the  beginning  of  occupancy,  *F. 

TDlv  -  Dry-bulb  temperature  of  inlet  air  at  end  of  each  hour  of  occupancy, 

•F. 

TD2  -  Initial  shelter  dry-bulb  temperature,  *F. 

T8S?K(l)  *  Initial  shelter  affective  temperature,  *F. 

TO  -  Initial  temperature  of  each  boundary  n?J*l  point,  *F. 

-  Total  radiation  trecssKtanre  of  boundary  y 

TZQltJ.j }  -  Teapera  tyre  of  outermost  nodal  paint  when  considering:  a  soil- 

backed  b-!\>ndary  j,  *F. 

ftvsi’KAfi  Mrt’i  r.(w<,uf.*J 
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TWB(l)  -Initial  shelter  wet-bulb  temperature,  *F. 


W  -Equivalent  interior  width  of  shelter,  feet. 

HCOIL  -Humidity  ratio  nt  saturation  corresponding  to  TCOIL, 

lb  noisture/lb  dry  air.  Suggested  value  is  0.00921  for 
TCOIL  =  55°F. 

Xllj)  -Inclination  angle  of  boundary  j,  measured  from  the  vertical, 
degrees.  Values  between  0*  and  $?0*  are  allowed.  For  a 
vertical  boundary,  Xl(j)  *  0*. 

Limitations  on  Input  Values 

Since  relations  exist  between  certain  inputs,  changing  one  necessitates 
adjusting  related  values.  These  related  input  values  are: 

1.  If  the  time  increment  DT  Is  varied,  adjust  IDATA(5),  since  IDATA(5)  a 
1  hour/DT,  rounded  "to  -Re  nearest  integer. 

2.  For  the  NW  boundaries,  the  boundary  properties  of  AREA(j),  ASOTL(j), 
AZ(j),  Kl( j),  K0( j),  LAYKAX(j),  Kt-!AX(J )  must  correspond. 

3.  If  the  number  of  layers,  LAYKAX(j),  is  changed,  then  layer  properties 
must  be  added  or  deleted:  CON,  CP,  DP,  IK,  ROW. 

4.  Equivalent  dimensions  CG,  HH,  W  must  correspond  to  the  ARFA(j)  values. 

5.  The  hour  limits  on  the  interval  for  which  average  temperatures  arc* 
desired,  IKT(rt.)>  must  not  exceed  the  length  of  occupancy,  HRS. 

6.  The  values  of  £(PP(,; }/DX( j ) )  ♦  1  must  equal  KJ'AX  for  each  boundary. 

i 

V.  Initial  values  of  dry-bulb,  vet-bulb,  and  effective  temperature  and 
relative  humidity  must  correspond:  TP-?,  TVR(l),  R?. 

8.  KSOU.  rust  c«. rresponA  to  TCOIL,  for  the  air  conditioning  unit. 


It  is  obvious  that  the  sire  of  the  data  deck  is  quite  large;  it  is  of 
ut-m-jft  imp r lance  -hat  the  order  of  the  data  dock,  as  indicated  in  Figure  3-2, 
to  »i*  -5.  pin*-?  av  t  ■c;d,hs  of  the  variables  are  also 


»'  untri.  A'*tri  *ah  *>;  r- micm  uiv^oh 


+ 


■•flL 

^jfT^yiTyJ 

Ifipi 

111 

I DATA  (all  13  elements) 


ll«ll«tdltMlltilllhlth«IMIlfcllb|g|||tatlblllllll||||||||||||||»ltl||ill|||lttllll 

in  hi  linn  mitt  limit  mu  tn  in  in  mi  ii  tit  tin  n  i  tint  i  ini  t  ii  ii  it  n  mm 


psat(i) 


(71  CARDS) 


ItltMdtMIMIIIIMmiMMUMIMMMIMliMMIIIIKtlimilimillilMllimil 

mi  in  iii  iiiiiiMiinii  iiitiinini  i  ii  mi  ii  i  ii  i  iii  nit  ii  im  it  t  ii  i  iii  1 1  in  ii  it 


KiTIH( I) 


(365  CARDS) 


Illicit  ItttHtlltltHtlllltttHIIHIttllttltlltlHItHIttltHMttttHllllt  MMII 
■  1  1  1  vi  ft  t  ««  • *  *»«  *•  »  ••»»*i»»iiiiii»aiB  a  »*  11  a»*  »»  »»»  »  •»• 

ttiiMniMiiniiiiitniiiiiitiiiiiiMiiiiiiiiiMiiitiiiiiiiiiiiiiiiiiiiiiiiiiii 


*  Bote:  The  #e$ucnce  of  input  card#  5-6- 7-8  i#  repeated  for  each  boundary 

in  succession;  i.e.,  for  J  fro«  1  to  KW. 

♦  Bote:  Input  card  6  is  included  in  the  r«ws<*«cc  for  boundary  J  only  if 

ASOrL(j)  <  0.0. 


Figure  3-?  Inp-is*.  r®lo  fcira* 
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/ 


DT 


12 


11 


V™ 


•  lltll  JtlMMIMMIMMIMMIMIMItMIMMlMMIMMIMMIMMIMIMMIIIMIttl 

mi  it  unit  hi  mm  iiiiniiin  t  mi  ii  tin  it  mi  it  iii!  iff  ttiiiiiiitniiiiiinn 


•  it 


ml 


^niiiiiiiMiiiiMiirmimmiiiMiiiMiMiiiiiiiiMiimiiMiMiiiiiiii 

i  •  i  •  a  aaaa»aaiiaaannnM»niia!t»»aH»aa»aa'a««aaaa«aa««BaM»aa»aaaii«iiaaiiaaa»BBaaa»aaM 


nniiiinint  until  t  mi  in  nniiiinit  nil  nun  in  nnni  nun nntinnin 


V- 


TS0IL(J) 


10 


mm(i  111111111111111111111111111111111111111111111111111111111111111111111111111 

>•!• ia«aaMaaaaaa»nnHaaaaaaiaaMaa»aaaaM«Maa««MM»ttuuaa*aaaauaHeai>aaaaaaaaanaita 

in  nit  t  nun  in  mi  nnn  nin  in  nn  in  mm  in  t  minn  nun  i  mi  ninnn 


CFMN 


>4 

M 

o 

8 


WCCJl 


BF 


F2 


g 


R?i 


im  inti  tjg  it  ni*Mm|Mai|omii!ii|MtJi*ao|iiti{n  iijnii  in  in  mi  in  miimiim 

ii><i<iiMiaalHaaaaaar«nMiiiapaaaiaas4aE*WBaac 


Is 


"N 


•  ■  RniURU  ft  *  H 


V, 


nnnnnn  tiinnnnn  in  mini  i  nni  tin  iniinn  nnn  iii  nnni  nni  inn 


DP  JJK  cm  $ai  CP 
(K,J)|(K,j)l(K,J)|(k,j)|(K,j) 


** 


inn  tfi  iim  Ji  nil  ill  nil  tji  nil  ilniiiiiiiiiiiiitiiiiiitiiiiiiiiiiiitiiiiiiiiiiiiii 

il  iiiaiiiaiaiXM  ah»iaaiaiaia»ibaa»oaR«>«i!«aoa«a>iiii»iiaaaMaaa>i>aa«i  aniMiata 

nntnmnnmnmmmmnmmmmtmmmmnmnmminmm 


/az(j) 


“N 


4-r 


•  iiiiitjnnnininnnnnnmninnnnnnnnininniinininninnnn 

■  mu  ■liMiMaiaMiaimaiia.iiiiiiiliiiaoiiaaaiiaaaiiuaiivaaaiacaiiiii'aaii'iaKMi'ia 
III  III  HI  III  III  III  III  III  I  III  III  I  I  III  III  III  111  I  III  III  III  III  III  III  I  III  III  I  I  III  I  III 


V 


*  See  Rote  on  preceding  page. 

44  Koto:  Input  card  7  is  included  in  the  sequence  only  if  AS0IL(J)  <  0.0 
and  XI(J)  <  90*0. 

«*  Rote:  Input  card  6  is  repeated  for  each  layer  of  boundary  J,  i,e.,  for 
K  free  i  to  IAYKAX?J)  for  boundary  J. 

:  Rote:  Input  card  10  is  included  only  for  those  boundaries  for  vhich 
AS0I1(J)  >  0.0. 

*  Kote:  Input  card  11  is  included  once  if  IIWTA(l)  0  .  If  IDATA(l)  -  1  , 

card  11  wist  be  repeated  for  each  nodU  point  of  each  boundary,  in 
succession;  i.e.,  for  K  lYow  l  to  LAYI4AX(j)  for  each  J  froo  1  to  W. 


Figure  3*?  (Confd)  Input  Data  Forrat 


IIHMIllMHIMMIIIMliniMIMIMMIIMIMIMIMMiMMIMMIIIIHIIIIIIIHH 
■  1 1  •  i  •  >  •  tuaaaaaaaa  aaaaaaaas Baa  t s  *«• Bn* aa *•« ««i»» aaaa a  aaaa 
tiiiniiniiiimnintiituiiMmiMmmimtttiinnnMittmiiiiiniinii 


IMIlliMMMMIIMIISMMMIIIIIISMntCIIIIMtIOMtfMMMMMIMHIIIIMIMM 
i  >  I  •  I  •  1 1  •  aa  aaaaaaaaaa  anana  it  ana  aaaa aaa  ikiiiomi  r  aaa Baa Baa  BMaaa.iJ  aauaaa  a  a  a  aaa  aaaa 

liiiiiiiimmiiiiiiiiiiiiiiiiiiiiiiiimiiimiiiiiiiii'miiiiiiiiiiiiiiiiiii 


§  3  § 
P  i  e  i  s 


i  >  i  <ii  •  »  *!«•?»  a  ota  aaaa  aoaBnt.B3BBBBi*auBBBBaaa«H*«iH«»*iB»auBBB»BMaBBOBBBWB»B»«o»Ba>t  ana 

iimiiiiiiiiiiiiitiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiimiiiiiiii]iiii!iiiiiiiiiiii 


CFKIN  QLTTE  QEQUXP 


IIIIMMIIMMMIMtlllHMttllMllltd 
l  «  « til  » •  I  iiaaalaaaaaaiagtiaisa»i«i»»»U) 


1 1  •  <  dl  t  •  ■  t  •  a  a  wa  a  a  a  aaa  ana?  a  a  a  bur  aaaa  aaa  aaWaaaaaaa  aaaa  ttHaivvattaaaaaaaaaaai  aa a  aitaaa 
1 1 1 1 1 1 1 11 1 1  It  II 1 1 1 1 II 1 1  ’  1 1 1  H  1 1 1  Ml  1 1  fll  1 1 1 II II 1 1 1 1 1  M  1 1 II  111  H  1 1 1 1  H  1 1 1 1 1 II II 1 1 


Miliii3  ta 


IlSlIlieiggHMIKtMlinitlKIIIMMIMIIMIMIMIMIIIMimilllMIIIIIMIIIIIIItl 

1 1  I  til  •  I  ill  i  i  ofa.a  J.aarrtoaaaiar . at  o»  aaaa  u  aa«  a.  .  b«b»  auaaaa  bbbbbi- a  nr  aaa  a  •  a  n  a  a  a  ri  ana 

II  I  1 1  I  II  I  I  I  I  I  I  I  II  I  I  I  I  I  I  I  I  II  1  1  I  I  1  I  I  I  I  I  1  II  I  I  II  I  I  I  I  I  I  I  I  I  I  I  II  I  II  I  I  I  I  II  I  I  II  I  I  I  I  1  It  I  I  I 


INT(I) 


M  *  •  *  •  •  •  II*  9  «»»»  ft  It  »•*»»»«*  <•«»••«««»»  VII  *«*«*  «••»■»%» 

n  1 1 1 1 1 1 1  m  1 1 1 1 1 1 1 1 1 1 1 1 1  n  i '  1 1 1 1 1 1 1 1  n  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1  n  t  n  1 1 1  m  1 1 1 1 1 1 1 1  n  i 


*  Note:  Input  card  13  it  Included  only  If  KOI NT  >  C,  and  is  repeated  for 
each  lover  and  upper  hour  nurfccr  for  each  interval  in  succession; 
i,«.,  for  I  fret*  1  to  2  s  (KOIKT). 

♦♦  Rote:  Input  card  15  i*  included  only  if  IDATA(12)  <  0,  and  it  repeated 
for  each  hour  of  occupancy  after  the  first; Ce.*  for  hour  nunbors 
frea  2  to  HRS. 

♦a*  Note:  Input  card  16  is  included  only  if  IBATA(13)  <  0,  and  is  repeated 
for  each  hour  of  occupancy  afte~  the  first;  i.e.*  for  hour  mnbers 
frees  2  to  I3ts. 


Fifurc  3-2  (Con’.'d)  Input  paU  Forr* 


indicated  in  the  figure.  Since  data  may  be  placed  only  in  specific  fields 
on  each  data  card,  the  easiest  method  of  keeping  the  data  deck  ordered  is 
to  punch  an  appropriate  identifier  anyvhere  outside  of  the  specified  fields. 

Values  for  the  matrices  BQTIM  and  PSAT  utay  be  found  in  Tables  I  and  II. 

The  EQTIM  values  are  those  for  the  year  1958,  and  may  be  considered  those  for 
any  year  in  this  century.  Values  of  the  saturation  pressure  as  a  function  of 
dry-bulb  temperature  are  introduced  as  a  table  with  a  linear  interpolation 
scheme  included  in  the  program.  These  two  inputs  comprise  the  first  U36  cards 
of  a  standard  data  deck. 

If  the  program  is  used  to  analyze  an  actual  shelter  under  actual  weather 

conditions,  hourly  ambient  temperature  and  relative  humidity  values  are 

- 

required.  Soil  temperatures  may  be  considered  constant  over  a  two-week 
occupancy  period. 

Various  time  increments  (DT)  were  tested  by  the  authors  to  determine  a 
compromise  between  accuracy  of  computation  and  length  of  computation.  It 
has  been  found  that  a  time  increment  ot  ten  minutes  (O.loY  hr)  produces 
results  which  are,  on  the  average,  within  0.3*F  of  the  results  calculated 
using  ? -minute  increments.  The  computation  time  is  also  decreased  by 
appryx irately  see  Figure  3-3* 

It  should  be  noted  that  integer  data  values  must  be  right-hand  justified 
within  the’r  respective  fields;  if  this  is  not  the  case,  the  machine  will  place 
■/erne.*  in  *  he  field  to  the  right  of  the  right-most  non-rero  digit. 


r  W  I./  r  I'lAM  iTv-ArCH  r.lVi  V  ' 
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Figure  3*3  Shelter  Effective  Temperature  Computed  on  Three  Different 
Time  Increments 


SECTION  4 


HOW  TO  USE  SHSP 

4.1  The  first  step  in  the  analysis  of  a  shelter  is  to  describe  the  shelter 
with  all  of  its  necessary  characteristics.  Answer  these  questions  on  a 
form  sheet  such  as  that  of  Figure  4-1. 

1.  Where  is  the  shelter  located?  (ALT,  ALAT,  ALONG ). 

2.  What  is  the  shelter  geometry?  How  many  separate  "walls"? 

(NW,  cg,  m,  W). 

3.  How  many  occupants  will  be  in  the  shelter  and  wi.il  this  number 
change  during  the  period  of  occupancy?  (P). 

4.  What  is  the  total  ventilation  airflow  entering  shelter  and  will 
it  change  with  time?  (Fl). 

5.  What  is  the  lighting  and  equipment  load  and  is  it  constant? 

(QLITE,  QJSQUIP) . 

6.  How  long  is  the  period  of  occupancy  and  when  does  it  start 
(day  of  the  year  and  time  of  day)?  (D#  HKS,  LST(l)). 

7.  For  each  boundary,  what  are  the  values  of  the: 

A.  Area?  (AREA). 

B.  Humber  of  layers  and  nurber  of  nodal  points?  (LAYHAX,  NKAa). 

C.  Interior  and  estcrior  film  heat  transfer  coefficients? 

(HI,  HO). 

P.  Control  number  for  the  wall  broking?  (AS01L). 

B.  Solar  radiation  transmittance?  (TRAM?)- 
F.  Inclination  and  azimuth  angles?  (XI,  A?.). 

6.  For  es~h  layer  of  each  boundary,  what  are  the  vaiu*:*  of  the: 

A.  Thickness?  (PP). 

B.  Kcdal  Pint  Spare?  (I*). 


BFKI  UAI  AMCIIICAM  Ri.T.rfttiSH  O iV;r:»O.M 


EQTIM1365) 


o::rjf;nAL  n*.v*»  Ari;.u  cjivi  j'o:: 


*1-2 


Boundory  AREA  LAYMAX  NMAX  HI  HO  ASOIL  TRANS 


fsi 


M 


n 


u 


n 


□ 


□ 


U 


u  u  u 


•J  *1  I  fl  I 


a  • 


1 

Li 


00  0 
0  0  0  0  0  0'0  0  0  0 


Z  "  N  10 


1 

0  U 


GGGGGGQOQGGQGGQGGQGQ 


0 


GGGQQGGGQGGGQGGQGGGG 

GGGQQOGGGGGGGGGGQQGQ 


i 

a 

i 

a 

• 

• 

I 

• 

• 

• 

a 

If 

• 

*j  r 

•i 

a 

a 

a 

10 

J 

• 

’ 

j 

: 

j 

11 

_»  L 

■ 

j  i— 

* 

I L 

• 

J  _ 

3 

• 

- 

3  li 

rriDtDf^OOOlO  —  PJfO%flf)(X)f^COo>0 

—  —  —  —  —  —  —  —  —  (\j 


*T 

<D 

i 

lO 


rjr.rjtrnAt  /.to:  t.iC(>rj  ri'.  ,t.  mi.:h 


4-3 


Figure  4-1  (Cont'd)  Form  Sheet  for  Organizing  Input  Data 
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Figure  4-i  (Cont’d)  Fora  Sheet  for  Organizing  Input  Data 


C.  Thermal  Conductivity?  tC0N). 

D.  Density?  (ROW). 

E.  Specific  Heat?  (CP). 

9.  If  air  conditioning  is  involved,  what  air  flow  is  going  through  the 
cooler  (CFKIN);  what  is  the  effective  surface  coil  temperature  (TCOIL) 
and  the  humidity  ratio  corresponding  to  this  dewpoint  temperature 
(WCOIL);  what  is  the  bypass  factor  (BP)?  What  is  the  amount  of 
outside  air  entering  the  system  (Ff*)? 

10.  Whet  are  the  initial  conditions  irsioe  the  shelter?  (TD2,  R2,  TWB(i), 
TEB’P(l).  What  are  the  initial  wall  temperatures?  (TO). 

11.  What  are  the  soil  temperatures  around  tr.e  shelter?  (TSOIL). 

12.  ?or  computational- purposes: 

A.  Hew  many  intervals  are  there  for  which  average  temperatures 
are  required,  and  vhat  are  their  hour  limits?  (NOHiT,  IhT), 

3.  What  tine  increment  should  shelter  temperatures  be  calculated 
on?  (W). 

13.  What  are  the  hourly  values  of  inlet  temperature  and  relative  humidity 

duri  ccupancy?  (TA2,  TD12,  R2).  What  are  the  values  of  o-cupar.cy 
variables  (P,  Fi,  ^KIK,  Q*$UJP;? 

After  the  dale,  is  determined,  options  in  the  program  must  be  decided 
upon.  Will  the  analysis  include: 

1.  Initially  isothermal  boundaries?  (IDATA(l)). 

?.  Adlabetic  or  non •c.d’.atatic  bour.darie??  (IDATA(6)}. 

3.  Air  conditioning!  (IDV.A{7)). 

5*.  Printout  of  th«’  tVclc  of  hourly  value*?  fPATA(ll)). 


5.  Constant  or  time-varying  ftccupency  variables?  (IBATA(12)). 

6.  Constant  or  time-varying  inlet  conditions?  (IDATA(13)). 

4.3  Once  the  data  has  been  determined  and  the  options  have  been  chosen, 
initial  conditions  in  the  shelter  system  must  be  set.  Four  quantities  must 
be  giver,  to  specify  the  shelter  atmosphere:  dry-bulb,  wet-bulb  and  effective 
temperatures,  end  the  relative  humidity.  It  is  mandatory  that  these  values 
correspond  to  each  other. 

Initial  well  temperatures  snould  approximate  the  initial  shelter  dry- 
bulb  temperature  within  a  few  degrees.  While  stability  of  solution  is  not 
u  problem  in  the  code,  initial  conditions  should  be  fairly  accurate  so  that 
the  number  o?  calculations  required  to  eliminate  the  errors  introduced  by 
the  incorrect  estimate  of  the  shelter  and  boundary  temperatures  is  minimized. 

4.4  After  all  values  necessary  for  an  onalysis  have  been  determined,  the 
code  user  is  ready  to  create  the  data  deck.  The  format  of  each  data  card 
and  the  order  of  the  cards  is  indicated  in  Figure  3-2-  Of  course,  the  order 
of  the  deck  is  of  major  importance  and  an  identification  on  each  data  card 
will  i.eip  keep  it  in  its  correct  position. 

4.5  SHKP  code  is  written  in  FORTRAN  IV  specifically  for  the  CDO  3600 
computer  system  and,  therefore,  utilizes  certain  library  functions.  These 

include: 

Form  Definition 

ABPr’(x)  Absolute  value  of  X 

ATAIi(X)  Aret.i\gc.ni  X  radian* 

COS(X)  Cosine  X  radians 

EXP(X),  KXFv(X)  e  to  the  X--  power 

Sllf(X)  Sine  X  radian? 

Ssi5T(x}  Square  root  of  X 

TA3i>‘(x)  rancor,*.  X  radians 

ertifi’M.  Afc.rrMRr.M  r.Arscw  nrv:  • 
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*  • 

If  the  code  is  re-written  for  another  computer  system,  these  library 
functions,  or  their  equivalents,  must  be  supplied. 

4.6  If  the  values  of  certain  input  data  and  certain  computed  variables  exceed 
specified  ranges,  error  messages  will  appear  in  the  output  as  per  Table  III, 
p.  2-l6. 

The  SHEP  code  requires  two  memory  banks  (o^  locations)  for  the  code  and 
the  accompanying  matrices;  in  fact,  the  AINV  matrix,  which  contains  the  inverted 
matrices  of  the  temperature  coefficients  in  the  boundary  heat-transfer  equations, 
fills  one  bank  (76,400  out  of  77,777  (octal)  words).  Thus,  a  BANK  statement 
is  another  requirement.  There  are  only  two  computer  systems  ether  than  the 
3600  that  have  this  amount  of  storage  known  to  the  authors  at  the  tine  of  this 
writing;  they  are  the  CDC  6600  and  the  UNIVAC  1108.  Although  several  smaller 
computers  can  be  adapted  and  enlarged,  this  requires  special  facilities. 
Therefore,  in  terms  of  storage  requirenents,  the  SHEP  code  can  be  run  on  ar.y 
of  these  three  computing  systems. 
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SECTION  5 


INTERPRETATION  OF  OUTPUT 

The  output  of  the  SHEP  program  is  a  tabulation  of  the  instantaneous  values 
of  the  thermal  loads  affecting  the  shelter  and  the  shelter  psychrometric  con¬ 
dition  for  the  last  increment  of  each  hour  cf  occupancy,  see  Figure  5-1 •  It 
must  be  remembered  that  these  values  are  instantaneous  and,  while  they  reflect 
trends  in  the  shelter  system,  interpolation  between  hourly  values  may  be  invalid. 
Also,  hourly  thermal  loads  determine  the  resultant  shelter  temperature,  not  vice 
versa.  The  3600  computer  system  printer  includes  a  procedure  for  rounding  output 
to  the  format  specified.  Thus,  the  thermal  loads  are  rounded  in  value  to  the 
nearest  integer. 

The  air  conditioning,  boundary  and  condensation  heat  loads,  when  positive, 
represent  energy  lost  from  the  shelter  system;  ell  other  positive  loads  are 
energy  additions.  In  the  SHEP  code,  the  assumptions  associated  with  the  cal¬ 
culations  of  the  air  conditioning,  moisture  condensation,  and  transmitted  solar 
radiation  loads  general] y  lead  to  conservative  estimates  of  these  loads. 

Maximum  and  minimum  shelter  temperatures  and  the  hour*  of  their  occurrence 
are  indicated.  Average  temperatures  are  computed  for  the  entire  occupancy 
period,  for  each  fu.il  day  of  occupancy,  and  for  specific  intervais  determined 
by  the  code  usar. 
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SECTION  .6 
SAMPLE  PROBLEM 

To  illustrate  the  use  of  the  SHEP  code,  a  shelter  defined  by  the  values 
listed  in  Table  IV  is  considered.  A  sketch  of  the  sample  shelter  is  shown 
in  Figure  601;  the  rectangular  shelter  area  is  in  the  basement  of  a  multi¬ 
story  building.  It  is  defined  by  eight  opaque  boundaries,  four  of  which  are 
multi-layered;  three  are  soil-oacked,  two  boundaries  are  exposed  to  ambient 
air,  and  the  remaining  three  are  exposed  to  interior  spaces.  See  Figures 
2-1  and  2-2.  Occupancy  is  limited  to  1000  persons  for  2k  hours  at  a  constant 
ventilation  rate  of  15  cfm/occupant.  The  air  conditioning  system  is  assumed 
to  be  operating  during  occupancy  with  a  jnstant  airflow  rate  of  8  cfm-occupant 
and  the  shelter  lighting  and  equipment  add  a  constant  thermal  load  of  11,250 
Btu/hr.  The  weather  affecting  the  shelter  is  that  which  occurred  a*  Montgomery, 
Alabama  during  1  August,  19^3. 

The  output  of  this  shelter  analysis  is  shown  ir  Figure  501.  It  can  be 
seen  that  the  shelter  temperatures  vary  in  a  diurnal  cycle,  as  would  be 
expected;  they  are  lowest  in  the  early  norming  and  arc  highest  in  the  late 
afternoon.  Air  conditioning  also  can  be  seen  to  vary  directly  with  the 
shelter  temperature.  In  the  cooler  hours  of  the  day,  i.e.,  from  10  PM  to 
9  AM,  the  shelter  boundaries  are  warmer  than  the  shelter  atmosphere,  and 
energ)  is  transferred  from  the  boundaries  to  the  shelter  air,  as  is  indicated 
by  the  negative  signs  on  the  boundary  heat  losses  for  these  hours.  During 
the  other  hours  of  the  day,  the  shelter  is  warmer  than  the  boundaries  and 
energy  leaves  the  shelter. 

If  the  value  c:‘  condensation  heat  loss  is  negative,  this  cor  resends  to 
evapora*. '  -s  rf  the  c.r.dersate  ;r.  the  vails,  and  -an  addition  of  energy  to 
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TABLE  IV 

SAMPLE  INPUT  DATA 
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TABLE  IV  (Cont'd) 
SA.'<TPLE  INFJT  LATA 
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Flcuro  6-1  Sketch  of  Shelter  Area  Used  In  Sample  Problem 


the  shelter  temperature.  Positive  value  of  the  condensation  heat  load  indicates 
moisture  is  condensing  on  the  shelter  boundaries  and  energy  is  being  removed 
from  the  shelter  air.  There  is  a  check  in  the  code  to  prevent  the  amount  of 
re-evaporated  moisture  from  exceeding  the  amount  previously  condensed. 

Metabolic  loads  are  predicted,  using  the  relations  developed  by  Houghten^, 
based  on  both  the  s'.elter  dry -bulb  and  effective  temperatures.  Since  these  vary 
by  only  a  few  degrees  during  occupancy,  the  total  metabolic  load  for  each  hour 
is  approximately  constant. 
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APPENDIX  A 


KSTZRMtHATIOtt  OP  SHElffBP  ENVIRONMENT 

Tine  study  of  the  Interrelationships  among  ventilation  rate,  thermal 
load,  and  boundary  losses  in  the  shelter  is  based  on  the  assumption  that 
the  shelter  atmosphere  may  be  represented  by  one  psychrometric  condition. 

A  shelter  model  that  will  take  into  account  the  thermal  effects  of  solar 
radiation,  air-conditioning,  moisture  condensation  and  re- evaporation, 
and  boundary  losses  is  developed  according  to  the  following  assumptions: 

1)  the  air  within  the  shelter  is  completely  and  Instantaneously  mixed; 

2)  the  film  heat  transfer  coefficients  are  constant  for  any  one 
boundary  surface; 

3)  the  radiative  energy  transfer  within  the  shelter  can  be  neglected; 

4)  the  condition  of  the  air  exhausted  from  tre  shelter  is  the  con¬ 
dition  of  the  shelter  atmosphere; 

5)  the  thermal  and  physical  properties  of  the  structural  materials 
end  of  air  are  not  temperature-dependent; 

6)  the  incident  solar  radiation  i *  absorbed  on  the  outer  surface  of  the 
boundaries  <*M  appears  as  conducted  energy  or  Is  transmitted  Into 
the  shelter  and  Is  considered  \r»  Instantaneous  load  with  the  lighting 
end  equipment  loads;  and 

7)  the  therm 1  loads  and  the  paychrometrlc  states  of  the  inlet  and 
exhaust  air  are  constant  over  short  time  intervals. 

The  shelter  model  is  therefore  a  volume  enclosed  by  boundaries  into 
which  sensible  and  latent  heat  loads  and  ventilatinc  Air  are  introduced 
and  free  vhich  eir  Is  exhausted  and  energy  is  lost,  tee  figure  A-l.  The 
gcvcrniiy;  equations  of  the  enclosed  volume  are  derived  free  the  conservation 
of  moss  and  the  conservation  of  energy. 
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Metobotic 
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►  A-i  Shelter  Tt<er&  Transfer 
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&  vv-^vv  ■ 0  <i> 

where  H  *  enthalpy  of  the  exchanged  air,  Btu 
Hg  ■  enthalpy  of  the  shelter,  Btu 
Qjl  «  metabolic  energy  of  occupants,  Btu 
^  ■  energy  due  to  shelter  lighting,  Btu 

«  energy  due  to  other  shelter  equipment,  Btu 
Qy  *  energy  conducted  through  shelter  boundaries,  Btu 
Q  *  energy  removed  from  she leer  system  by  air  conditioner,  Btu 

Av 

Integrating  equation  (l)  ovor  the  time  increment  ft  yields 

•  V*  v  vv  vv  ■ c  (2> 

where  the  subscripts -1  and  2  refer  to  tbe-V  '  'ig  and  the  end  of  the 


time  increment,  respectively.  Althov 
in  the  sense  that  results  ur’ lg 
changing  conditions,  by  assvr.  • 
psychrometric  states 


considered  < 


lysis  is  time-varying 
increments  represent 
crmal  load  values  and  the 
in  this  equation  are  all 
.rsereforc,  there  Is  no 


beat-storeg  terr 
Conservation 
water  vapor  in  the  sh«l< 


-i  8.  V  Jp 


d  C  -  0. 


ooth  the  dry  air  and  the 


Vhen  it  is  applied  to  the  vapor; 


&  ^  ^,8'  %  'VcOKIl  ‘  ° 


vbere  Ny  ■  ***•  of  water  vapor  in  the  exchanged  air,  lb. 

•  mass  of  water  vapor  introduced  by  the  shelter 
**"  occupants*  lb. 

My  #  •  mass  of  voter  sapor  in  t»>*  shelter  air,  lb. 

My  cogs  *  mass  of  vator  vapor  removed  by  condensation,  lb. 

Integral?  n  ylcld> 

^*V,l  *  '  *Sf,S,2^  4  *V,M  '  N,C0i 3>  "  0 
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The  energy  balance  for  dry  air  in  the  eysteo  ie  given  by 


&  <M.>  *  &  ■  0  <*> 


where 

* 

V.  - 

oass  of  d-y  ail  in  the  exchanged  air,  lb. 

mass  of  dry  air  in  tl*  shelter  air,  lb. 

or 

<Vi 

•  Vo’  *  (1W  ‘  V.,a>  •  0 

(6) 

let 

6  Q  « 

WWW V 

(7) 

V« 

*  Vi  *  Vs,i  *  Vm  "  Vcoro 

(8) 

V» 

*  *  ,  +  ><  , 

a,l  .  a,e,3. 

(9) 

then 

6q  ■ 

h*\t 

(10) 

Vo 

"  Nfa  *  \s,8 

(U) 

M*>o 

*  \,2  *  \j9,2 

(12) 

Nov 

*  Ma,8,2  ha,M  *  *¥,8,2*7, 3, 2 

(13) 

and 

*  • 

*a,2  h -A  *  V*  *v,2 

(1*0 

vhere  the  lover  cma  h's  refer  to  specific  enthalpies  (Btu/lb). 


Since  the  ventilating  air  lccvec  rt  the  shelter  condition, 

\ ,2  " 
and 

Va  *  V$,£ 

Substituting  this  value  for  the  i ncrascntal  energy  6 

**  *  \,0  h»,«,8  +  *V,0  N,!,8 


(15) 


(16) 


(17) 
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lev  specific  enthalpies  of  dry  air  >  A  water  vapor  a>y  be  given  as 
functions  of  dry-bulb  tempers.tuie: 


0.24  Tdb  (18) 

»  for  32.2*?  <  £  150 *F 

1061,0  ♦  0.444  (Tdb)  (19) 


Combining  equations  (17)  through  (19)  and  solving  for  the  dry-bulb 
temperature  leads  to 


6ft  -  1061  Hyfo 

0.444  ♦  0.24  \  0 


(20) 


This  represents  the  shelter  dry-bulb  temperature  at  the  end  of  the 
tine  increment  6r> 


Expressions  for  &  and  H^Q  can  be  found,  since 

M..a  *  60  V.,a  <T 

(21) 

Va  ■  60  V».«  4T 

(22) 

t.i,a  * 

(23) 

“4  V»,8  *  V',».S1! 

(24) 

vh*ra  V  ■  ahalttr  voXwm,  ft4 

p  _  •  density  of  dry  air  leaving  shelter,  lb/ft^ 

2  *  density  of  water  vapor  leaving  shelter,  lb/ft ^ 

p  a  »  density  of  dry  air  in  shelter,  lb/ft ^ 

***** 

.  .  »  density  of  water  vapor  ir.  shelter,  lb/ft^ 
and  ^  *  mass  flow  rate  of  exhaust  air,  lb/hr. 
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Ve  bm  assuaed  thit 

P|(2  *  (25) 

aRd  *v,2  *  Pv,«,2  <26> 

therefore,  we  obtain  the  expression: 

\0  “  *>v,2(6°  VT+V)  <*> 

•“4  M^0  -  2  (60  P2  6t ♦  V)  (28) 

• 

Bow,  when  these  equations  ere  solved  for  the  ventilation  rate  fg, 
i  .  Vo  *_Pa,2V  1  Vo  "  PV,2V 


P2  1 


rPa,2 


>v,2 


It  is  assumed  that  the  densities  of  dry  air  and  c f  water  vapor  aay  be 
detemlned  frees  the  pcrfset-fAS  law:  Therefore, 


V  * 


?b*Pv 

R?W1 


O 

where  *h  *  barowotric  pressure,  lb/ft 

2 

Py  •  partial  pressure  of  water  vapor,  lb/ft 
8y  *  fas  coos tent  for  water  vapor,  ft-lb/lb  Dole  -  *8 
and  •  fas  constant  for  dry-air,  ft- lb/lb  sole  •  *8 

with  •  (9.8  *  10*7)  A*  -  0.0799  A  +  2U6.8  (38) 

where  *A  a  altitude,  ft. 

Bow,  the  partial  pressure  of  watet  vapor  is  expressible  at 

V '  •  r  \  (S3) 

where  r  *  relative  hsaidUy  in  the  shelter 

end  Pg  «  saturation  prta«<rn  of  water  Vapor,  lb/ft* 
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Saturation  pressure  varies  with  dry-bulb  temperature  according  to 


Pg  -  5.132  exp  (0.0329  Tib) 


m 


Using  to  values 
R  «  5J.35  ft- lb/lb  mole  -  *R 

ct 

and  Rv,  «  85.71  ft-lb/lb  mole  -  *R 
leads  to 

P,  •  r  (5.132  exp  (0.0329  Tdb)  )  (35) 

py  *  7V1"  -/lisgsgj  <°-ws87  «*p  <0.0329  Tab)  (36) 

30 

and  pa  *  t  -~\3q.69  **  1‘64  *•  1Q-‘8  “  0.001U2A 
db 


To  aswirtie,  the  tuo  equation*  nec^rsarj  to  coelute  the  theitnr 
cond-t^on  are  (20)  and  (to).  The  quantities  in  these  equations  **y  v* 
e<stpu*ed  as  outlined  above.  Although  the  dry- bulb  temperature  and 
relative  lucidity  uniquely  determine  the  condition  of  the  shelter  air, 
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it  is  desired  to  express  this  condition  in  terms  of  dry-bulb  temperature 
and  vet-bulb  temperature. 


The  wet-bulb  temperature  car.  be  determined  using  the  Carrier  equation: 


,<pb  -  *.  >  <Tab  -  T»b> 


P  P  /  W  »  UU 

*V  "  3  *  '  2d00.0  -  1.3  T 


vb 


) 


m 


where 


'  •  2 
P  *  saturation  pressure  at  wet-bulb  temperature,  Ib/ft 


vb 


vet-bulb  temperature,  *F 


P_  may  be  evaluated  from  equation  (3*0  as 

O 

Pg‘  <**  5-132  exp  (0.0329  Tyb) 


(42) 


Py  is  obtained  from  equation  (39)*  Equation  (41)  is  then  a  transcendental 
equation  in  Tvb  and  may  be  solved  by  an  iteration  procedure.  After  both 
the  dry- bulb  and  vet-bulb  temperatures  are  known,  the  effective  tempera¬ 
ture  may  then  be  calculated  from 


ET 


107.5  Tdb  -  45.2  Tvb 

62.3  ♦  T 


for 


db 


vb 


45*F  £  Tdb  S  110 *F 
3°*F  S  lwb  £  100 *F 


(43) 
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APPENDIX  B 

DETERMINATION  OF  SOLAR  RADIATION  INTENSITY 

B-l  Determination  of  the  Sun's  Altitude  and  Azimuth  Angles 

Consider  an  earth-centered  right-handed  rectilinear  coordinate 
system  with  the  Y-axis  directed  toward  the  sun  and  with  the  X-Y  plane 
that  of  *he  equator,  see  Figure  B-l.  The  unit  normal  vectors  "i,  J*, 
and  k  are  parallel  to  the  X,  Y  and  Z  axes,  respectively.  The  rays  of 
the  sun  are  parallel  to  the  unit  vector  s.  The  unit  vector  n  at  a 
point  on  the  surface  of  the  earth  P  is  normal  to  the  position  vector 
of  P  and  oriented  toward  the  north  pole.  The  position  vector  is 
parallel  to  tne  unit  vector  r.  The  point  P  is  located  on  the  surface 
of  the  earth  by  the  hour  angle,  6  ,  and  the  latitude  angle,  ^  . 

Now  r  *  cos  sin  ©  i  -  cos  $  cos  ©  J  +  sin  1c  (l) 

*  =  cos  7}  j  +  «in  T)  £  (2) 

whera:  Tj  °  declination  angle  of  the  sun 

The  sun's  renith  angle, f  ,  is  given  by 
*  arc  cos  (  F  *  i*  ) 

w  arc  cos  (sin  $  sir.  T)  -  cos  $  cob  ©  cob  n  )  (3) 

The  sun's  altitude  angle,  *3  .  la 


P  *  f  -  if 


Thuc 


0  *  arc  sin  (*'n  sin  T)  -  cot  $  coc  ©  ccsr))  (5) 


Now  the  vector  n  )■ 

«*  -sin  $  sin  0  i  ♦  sin**  cos  ©  J  *  cos  #  k 

TuMrg  R  •  *,  since  *  •  r  ■  0,  gives 

-•  .  - 
n  •  s  *>  sin  ©  cos  •  cos  *  cos  *'  sin  tj 


(6) 


(7) 
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Also  n  •  s  ■  cos  T  cos  |3 

where:  T  ■  sun's  azimuth  angle. 


Therefore, 


T  *  arc  cos  [sec  /3(cos  sin  7)  +  sin  $  cos  0  cos  7))]  (9) 


Thus  the  sun's  altitude  and  azimuth  angles  are  determined  as  a  function 
of  latitude,  $  ,  solar  declination  angle,  rj  ,  and  time  angle,  9. 

B-2  Solar  Declination  Angle  and  Time  Angle 

The  solar  declination  angle,  tj  ,  is  given  as  a  function  of  the  day 
of  the  year,  see  Figure  B-2.  This  variation  was  approximated  by  the 
following  function,  namely, 


V  =  0.410  sin  [  IBS75  (  -^o)  1 

where:  D  *  day  of  the  year,  1  to  3^5. 


The  tire  angle,  9,  is  denoted  by 


&  {  Tsolar) 


where 


T  ,  *  solar  time,  hrs,  see  Section  2.2. 


B-3  A’-glc  of  Surface  to  Sun's  Rays 

Consider  a  tilted  plane  surface  at  an  angle,  |  ,  to  a  vertical 
surface  at  a  point  P  on  the  surface  of  the  earth,  see  Figure  B-3-  The 
unit  normal  to  the  tilted  surface,  is  *t  tho  angle  |  to  the  normal 
to  the  vertical  surface  which  makes  an  angle, T  -  a,  with  respect  to  the 
projection  of  the  spin's  rays  on  a  horizontal  surface  at  point  P.  If  a 
right-handed  rectilinear  coordinate  sysla*,  is  defined  at  point  P,  such 
that  the  X-axis  i'  parallel  to  the  direction  of  the  projection  of  the  sun's 
ray,  a,  on  the  horizontal  plane  and  the  Y-axis  is  parallel  to  the  vertical  plane, 
the  X-T.  plane  is  the  horizontal  plane,  I/:t  tho  unit 
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Figure  b-2  Variation  of  Declination  Angle,  r,  ,  vith  the  Day  of  the  Year 


Figure  E-3  Geometry  for  Ti 


vectors  ~i,  J  end  k  be  directed  along  tha  X,  Y  and  Z  axes,  respectively. 
Then  the  unit  normal  to  the  tilted  surface,  is 

n  ■  cos(T-a)  cos  5  T  ♦  sin  5  "j  +  sin(T-a;  cost  'k  (12) 

3  s  a  «  azimuth  of  the  normal  to  the  vertical  surfs  measured 
westward  from  north 

T  *  azimuth  of  sun  measured  westward  from  north 
with  T  »  cos  fit  +  sin  /3j* 

Then 

s  •  n  ■  cos  fl  cos(T-a)  cos  %  +  sin  fl  sin  4  (13) 

but  7  •  n  »  cos  a  (14) 

B 

where:  -  - 

a  =  angle  between  the  normal  to  the  tilted  surface  and 
the  sun's  rays 

Therefore 

a  ■  arc  cos  (cos  0  cos(T-a)  cos  £  +  sin£  sin  4  )  (15) 

B-4  Sunrise  and  Sunset  Values  of  Hour  Angle 

At  sunrise,  /5  *  0;  then  from  equation  (5)  have 
0  «  arc  sin  (sin  d>  sin  -  cos  $  cot  cos  rj  ) 
or  «  arc  cos  (tan  $  tan  rj  )  (16) 

where:  Qr  »  value  of  0  at  sunrise 

However,  the  relation 

*1  £  tan  $  tan  T)  <11 
must  hold. 
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How  the  calculation  of  solar  radiation  is  restricted  to  tho 
northern  hemisphere! 

0  «  <  J  (17) 


and,  from  Figure  B-2,  T)  varies  as 
-0.410  <  73  S  0.410 


(18) 


Thus,  to  satisfy  the  restriction  on  the  $  and  tj  values,  the  ranges 
must  be 


$  +  V  £  \ 

and  *  £  | 


(19) 

(20) 


Since  7]  is  a  function  of  the- day  of  the  year,  these  restrictions' 
indicated  that  there  is  a  maxim®  value  of  latitude,  ,  that  can  be 
considered  at  ?ach  day  of  the  year,  see  Figure  3-1. 

The  value  of  the  hour  angle  for  sunset,  ©#,  is 

0,  «  2D  -  9r  (21) 

Iv'-s .  for  0  <  0  <  0  there  can  be  sor.e  fona  of  solar  radiation  upon 
r  -  -  s 

a  surface  as  long  as  0  <;  £  <  ^  . 

B-5  Intensity  of  Solar  Radiation 

The  intensity  of  solar  radiation  that  will  be  incident  upon  a 
surface  is  given  by 


I 

where  ? 


-  I 


DIR 


*DIR 

*DfF 

*REF 


4  XDIF  +  W 


(») 


•  Intensity  of  direct  solar  radiation 

•  intensity  of  diffuse  solar  radiation 

•  intensity  of  solar  radiation  reflected  free*  the  ground 
and  surrounding  surfaces  onto  the  surface 
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The  Assumption  is  made  that  th*  reflected  solar  radiation  from  the 
surroundings  is  negligible,  1^  ■  0. 


Thus 


1  “  rDIR  +  IDIF 


(23) 


The  direct  solar  radiation  is  given  by 


DIR 
where :  I 


1DN  C0“  tt 


0  <  a  £  */2 

a  >  2 

2 


(24) 


*  intensity  of  direct  normal  solar  radiation \  i.e.,  the 
DW  direct  radiation  on  a  surface  normal  to  the  sun’s  rays. 

The  intensity  of  the  direct  normal  solar  radiation,  1^,  is  given 

as  a  function  of  the  sun's  altitude  angle,  £  ,  (from  Reference  5),  see 

Figure  B-4.  The  dvta  &r.  approximated  by  the  relationships 

-  - 

294  (1-e  ^  )  for  clear  air 

lnri  “1  “  n  P  (25) 

* 20C  (1-e  )  for  industrial  air 


for  p>  >  0. 


The  intensity  of  the  diffuse  solar  radiation  is  assumed  to  be  defined 
by 


IDIF  *  XDIF,V 


where : 


DIF,  V 


IDIF,H  " 


+  ^DIF,  H  "  IDIF,V^  ,in  * 

intensity  of  diffuse  solar  radiation  on  a 

vertical  surface 

intensity  of  diffuse  solar  radiation  on  a 
horlsontal  surface. 


(26) 


The  ratio  of  the  intensity  of  diffuse  solar  radiation  on  a  vertical 
surface  to  the  intensity  of  diffuse  solar  radiation  on  »  horlsontal 
surface  is  given  as  a  function  of  the  cosine  of  the  angle  between  the 
norm!  to  the  vertical  surface  and  the  sun's  rays  (flroo  Reference  7), 
cee  Figu-s  R-5. 
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Normal  Radiation  as  a  Function 


The  results  are  approximated  by 


0.44  for  cos  $  cos  (?-a)  <-0.4 

*  0.625  l  cos  P  cos  (T-a)]2  *  0.52?  cos  0  cos  (T-a)+  0.S5  (27) 


DIF,H 


9  for  -C.4  <  cos  P  cos  (T-a)  ) 

^0.35  (  cos  £  cos  (T-a)  ]  c  +  0.45  cos  P  cos  (T-a)  ♦  O.55 

for  cos  £  cos  (T-?)  2  0 


SlF  "  Iwp,h  1(1  *  Bin^f  U  *T-&)  ♦  »in  •;  3  ^8) 

The  ratio  of  the  total  intensity  of  solar  radiation,  1^,  on  a  horizontal 
surface,  to  the  Intensity  of  diffuse  solar  radiation  on  a  horizontal 
surface,  H  is  given  as  a  function  of  the  diiect  normal  solar 
radiation  intensity  ( from  "Ref erence  y',  see  figure  B-6,  The  relationship 
can  be  approximated  by 

^22 -  -  0.667  *  10‘6  (Ij^)3  -  1.5  x  10~4  (Im)2+  °*ojl83  ( Ij^j)  ♦  1  (29) 

DIF,H 

It  may  be  assumed  t'oat 

XTH  "  XD1R,H  +  XDIF,H 


13W  *in  ^  +  IBIF,K 


Therefore, 


I  * 

W,H 


•w  ,l0  f 


DIF,  M 

Conblring  equation  (31)  and  equation  (20)  gives 

I  sin  p  -*)  (1  ■»  sin  O  ♦  *in  t  ) _ 

0ir  0,667  x  10*®  (Ijjj)*  -  1  5  x  lO*4  (Iw>4  O.OI83 


vith  iota),  radiation  on  a  surface  given  by  equation  (??}- 
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APPENDIX  C 

BOUNDARY  ENERGY  TRANSFER  AND 
NODAL  TEMPERATURE  COEFFICIENTS 

Heat  transfer  by  conduction  through  the  shelter  boundaries  is 
determined  by  the  temperature  of  each  boundary  inner  surface.  To  calcu¬ 
late  these  temperatures,  the  temperature  profile  through  each  boundary 
is  determined  by  the  solution  of  a  set  of  implicit  finite-difference 
equations.  Since  each  boundary  nay  be  multi-layered  (up  to  5  layers), 
the  nodal  poirt  coefficient*  must  be  determined  individually  for  each 
layer . 

For  this  solution,  a  l*»yer  is  divided  into  slaks  6X  thick.  A 
temperature  nodal  point  is  assigned  to  the  node  of  each  slab;  for  the 
interior  slabs,  this  point  is  at  the  midpoint  of  each  slab,  and  for  the 
innermost  and  outermost  slabs,  there  is  another  point  at  the  surface  of 
the  slab.  If  the  boundary  under  analysis  is  multi-layered,  a  similar 
nodal  point  array  accompanies  each  layer. 

Nodal  point  temperatures  are  determined  by  the  appropriate  heat 
transfer  equation  at  each  nodal  point.  There  are  four  basic  heat 
balances  that  can  occur  at  a  node;  the  one  to  be  used  depends  on  the 
position  cf  the  node  in  question.  These  four  positions  sc  (l)  at  the 
Inner  mrf&cc  of  tha  boundary;  (2)  in  the  interior  « f  any  layer  of  the 
boundary;  (3)  at  an  interface  between  layers  of  the  boundary;  and  (b) 
at  the  exterior  surface  of  the  boundary,  sec  Figure  C-l. 

FcriUon  3  -  nt  the  inner  surface  of  a  boundary.  The  heat  transfer 

in  this  case  is  by  convection  with  the  shelter  air  and  by  conduction 
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Interior  Interface  Exterior 

Surface  Surface 

Figure  C-l  Various  Positions  of  Boundary  Nodal  Points 


with  the  next  tlab  in  the  layer.  Heat  transfer  due  to  moisture  conden¬ 
sation  may  also  be  present,  see  Appendix  D.  The  general  heat  balance  for 
deriving  the  various  nodal  point  temperature  equations  is 
energy  entering  slab  n 
frem  preceding  slab  (n- 


energy  leaving  slab  n 
to  succeeding  slab  'n+l) 


energy  3tored 
in  slab  n 


For  position  1,  the  equation  becomes: 

^condensation  +  ^convection  ’  ^conduction  *  ^stored  ^ 

or  C  ♦  h,  -  Tn)  -  (Tn  -  Twl)  =  (T„  -  T„')  (2) 

where  C  =  condensation  term,  ste  Appendix  D. 

—  *  2  * 

*  interior  film  heat-transfer  coefficient,  Btu/hr-ft  -°F 


T 

db,s  =  dry-bulb  temperature  of  shelter,  #F 
Tr  =  temperature  of  node  n,  *F 

k  *  thermal  conductivity  of  layer  material,  Btu/hr-ft  -*F 
fix*  thickness  of  slab  in  layer,  ft. 

*  temperature  of  node  (  n+l),  *F 
p  =  density  of  layer  material,  lb/ft^ 

0^  «>  specific  heat  of  layer  material,  Btu/lb 
fi T  *  time  increment,  hr. 


Tr  *  temperature  of  rode  n  in  preceding  time  increment,  *F. 

The  matrix  solution  of  simultaneous  equations  requires  that  the 

coefficients  of  the  nodal  temperatures  arc  dimension1 ess.  Therefore, 

the  transformed  energy  balance  for  the  inner  surface  n>l*i  point  is 

%  pc  fix  .  pefix  , 

<i  0  *  KJTS  *  Fsflr  !  Tr,  •  <i CS~>  Tn.J  *  FiyS?  Tn 

*  V.  *  K7  (3> 
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which  in  matrix  element  notation  is 


fcnnTn  +  &n,n+l  Tn+1  “  bn 

cr,  since  the  inner  surface  nodal  point  is  indexed  1, 
all  T1  V  a12  T2  "  bl 


Position  2  -  in  the  interior  of  any  layer  of  the  boundary.  In  this 
case,  conduction  into  and  out  of  slab  n  is  the  only  heat  transfer.  There¬ 


fore,  the  energy  balance  becomes 


h  <Vi  -  V  -h  -  W 


pc&x 


<Tn  -  Tn  ) 


where  T  ^  ■  temperature  of  node  (n  -  1),  *F 


This  becomes,  in  dimensionless  quantities, 

(dx)2p 


■Vi  +  ("Or  +  2)  Tn  * 


pc  («*r  . 

- rrr -  T 


which  is  of  the  form 


*n,n-l  Tn-1  +  an,n  Tn  +  4n,n+l  Tn+1 


Position  3  -  at  an  interface  betveen  layers  (a,  and  (m+l)  )  of  a  boundary. 
Again,  the  only  heat  transfer  is  conduction,  but  since  the  successive  layers 

will  have  different  thermal  properties,  the  energy  balance  is: 

,  ,  p  c  6 

K  .  kfc+l  ■  *>'“ 

<Vi  -  v  -  <*.  -  w  ■  r— " 


or,  in  diner.s’onless  form, 


k  ,6 


x  p  c  $x  )  p  ,ep  ,6*  6* 

B  .  •  p«  m  ^  ro*l  s>  « 


»  .  #,  .  »*1  B  .  »  i*  »  _  rB*l  8*  1 

•  VX  *  ^  ^ -rs^r 


♦  (-  1  t  .  -  f  - -*U - _ - *♦-*  N 

^vT  nU  ?kBrr  ) 


p  k  rr 

a 


ecscsAi  nMraicAN  rfmarch  o*vtsto*>< 
C-‘» 


which  is  of  the  form: 


a  ..  T  ,  +  a  T  +  a  T  ,  ®  o 

n,n-l  n-1  n,n  n  n,n+l  n+1  n 


Position  4  -  at  the  exterior  surface  of  the  boundary.  This  case  has 
several  variations,  depending  on  whether  the  external  surface  boundary 
is  exposed  to  either  soil  or  air. 

If  the  boundary  is  backed  by  soil,  the  only  mode  of  heat  transfer 
will  be  conduction.  The  only  difference  between  position  4  and  position  2 
is  that  the  unknown  temperatures  T^+^  in  equation  (7)  is  replaced  by  the 
known  temperature  of  the  soil,  T  and  this  term,  being  known,  becomes 
part  of  the  b  coefficient: 


Pc  (6x) 

Vl  * 


pcp(6x)‘ 

kor 


T  +  *T 
ln  +  soil 


If  the  backing  of  the  wall  is  ambient  air,  convection  enters  as  well 
as  solar  radiation.  The  energy  balance  is 


(T  .  -  T  )  ♦  q  ,  -  h  <Tn  -  Ta) 

Ox  n-1  n  ^solar  o 


pc  fix 

*  v\ 


rr<T„-Tn>  (13) 


where  ‘Igolar  *  solar  radiation  load  on  boundary,  Btu/hr. 

h  **  exterior  film  beat  transfer  coefficient,  Btu/ftr-ft  -*>' 
o 

T  *  dry-bulb  temperature  of  ambient  air,  *F 


In  dimensionless  fern, 


(  'hT^n.1  4  IbT*  *  1  +  OF  1  Tn 


Both  equations  (12)  and  { lU ).  are  of  the  forn 


pc  fix  , 


4  Ta  4  SnT  Tn 
o 


a  ..  T  ,  ♦  a  T  •  b 
n,n-l  n-1  n»n  n  n 
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C  0  ard  T  is  the 

If  the  «01  i»  bached  *>  lnterl°r  *  ’ 

of  the  shelter  temperature  -  -  — * 

Kter  .u  coefficients  of  -X  nodal  point.  U  .  —  - 

«  .  «e  solved  simultaneously  vein*  a  Cau.e- 

>*'“  evaluated,  the  e*»a  on  ^  tBraugh  a  boundary 

Jor<te„  method  of  matrix  Inversion.  Convent 

is  then  computed  as  ^ 


<ib  ■  hl  tTdh,e"  V 


,  Mrtture  cf  the  inner  surface  nod.)  point  and  the  total 
where  T}_  *  temperature  oi 

boundary  heat  loss  is 


(YT) 
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APPENDIX  D 


MOISTURE  CONDENSATION 

T ha  process  of  moisture  condensing  on,  or  re-evaporating  from,  a 

surface  involves  both  a  mass  and  an  energy  transfer.  The  analysis  used 

( 9) 

in  this  study  parallels  that  of  Jakob  . 

The  following  differential  equations nay  ue  set  up  for  moving  fluids, 
relating  thermal  conduction  and  mechanical  diffusior: 


dT  .  v  3t  3t  .  v  dr 
3t  vx  3x  Vy  by  vt  31 


and 


|e.v  & 

ot  x  ox  y  oy  z  z 


«  (^£ 

3x2 


3y2 


dZT) 
+  3?) 


*3 


where 


T  ■  absolute  temperature,  *F 
v  «  velocity,  ft/hr  g 

p  «>  partial  pressure  of  water  vapor,  lb, 'ft 

'» . 

a  ■  tnermal  diffusivity,  ft  “/hr 
6  «  mechanical  diffusivity,  ft2/hr 


0) 

(2) 


Newton’s  law  yields 


and  Fick’s  law  gives 


in  •  ^ 

*  " 


where  q" 

m 

rate  of  heat  flow  per  unit  area,  Btu/ft  -hr 

m 

g 

rate  of  mass  flow  per  unit  area,  lb/ft  -hr 

n 

m 

direction  of  heat  and  mass  flow,  ft 

k 

m 

thermal  conductivity  of  the  medium,  6tu/hr-ft-*F 

c 

m 

concentration,  lb/ft^ 

(3) 

(«0 
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where  Ax^  =  mechanical  boundary  layer  thickness,  ft. 


and 


VQ  *  humidity  ratio  of  shelter  environment,  lb ^/lb^ 

Wg  *  humidity  ratio  at  boundary  temperature,  lb ^/lb^ 


Parallel  expression  to  equations  (10)  and  (11)  are 

q  «  h  (T0  -  Tg)  (12) 
•«d  \  -  pt b  (WQ  -  W8)  (13) 
where  b  *  t  coefficient  of  mass  transfer,  which  must  be  evaluated. 


A  comparison  of  equation  (10)  with  (12),  and  (11)  with  (13)  shows 

h  *  (14) 


k 


and 


6 


(15) 


It  is  assumed  that  the  thermal  and  mechanical  boundary  layer  thickness  are 
very  nearly  equal: 


AXT  = 


■  Ax. 


Therefore: 


or 


h 

b 


k 

s 


a 

1 


where 


Therefore, 


neb  p  c  6 

P  P 

a  *  ' thermal  diffusivity  * 

b  .  (-i-) 

Vp  a 


(16) 

(17) 

(18) 


P  e 

•  P 


(19) 


How,  for  water  diffusing  into  *ir,  ^  «  0.917  (Ref.  9) 

since  a  *  0-070  ft?/hr  for  air 
and  d  «■  0.?6*  for  water  vapor 
at  0*C  and  ?6C  m  Kg. 
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Therefore 


b  *  77P  otfiT' 


and 


1-09  (jV> 

*  P 


4v  *  P.b  («„  "  W,> 


■  1-09  (rV)  (K  -  W  ) 

'a  p 


vhere 


0^  *  specific,  heat  of  air-water  vapor  mixture 


*  0.2U  Btu/lb-’F 


Since 


q.  *  ny\ 


where  X  *  heat  of  condensation,  Btu/lb 


q  *  (v  -  W  ) 

^condensation  pc  '  o  a7 
^a  p 


Nov,  since  W  *  0.622  — 


Pb  *  P 


where  W  *  humidity  ratio 

p*  *  partial  pressure  of  water  vapor, 
p^  »  barometric  pressure, 
we  may  approximate 


V  -  V.  =  0.622  (  ) 

°  ‘  VPb-P0 


o  b  *  i 

Therefore  the  heat  transfer  d‘;«  to  condensation  is  expressible  by: 


Xh  ,  Po  P« 


) 


(2C 


(23 


(22 


(23 

(2l( 


(25 


(26: 


and  condensation  is  determined  to  be  present  if  pQ  >  n 


A  '  Po  Pb  '  Pw 

This  is  the 

condensation  term  of  equations  (1)  and  (2)  of  Appendix  C.  Also,  equation 
(?6)  expresses  the  heat  transfer  due  to  re-evaporation  of  the  condensate, 
If  the  energy  and  mass  transfer  phenomena  are  attuned  to  fee  equivalent 
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in  the  two  processes.  In  the  cftse  of  evaporation,  p  >  p  and  the 

S  0 

Condensation  terra  i9  ce«ative  (c  <  °)-  The  vaiue  of  Pe  in  elation  (26) 
is  assumed  to  be  defined  by 

pa  -  3.25  x  10"7  Td  -  1.86  x  10"4  T  +  O.O863  (27) 

where  T  ■  the  average  of  the  temperature  of  the  boundary  surface 
in  question  and  the  shelter  dry-bulb  temperature,  *P. 
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where  *  radius  of  curvature  of  the  inner  surface  of  the 

boundary 

■  radius  of  curvature  to  the  nth  nodal  point 
See  Figure  E-l. 

This  is  equivalent  to 


where  A 

E 

and 

0  represents  the  difference  in  calculating  heat  transfer  using  the 

two  analyses ;  it  represents  the  percent  departure  of  the  one-dimensional 

analysis  from  the  three-dimensional  analysis.  Values  of  0  as  a  function 
r 

of  —  are  shown  in  Table  V. 
rl 

Using  a  similar  analysis  for  cylindrical  coordinates,  the  expression 
for  0  becomes  r 

0  -  $  (r~~)  ln^f)  (7) 

-s-  i  1 

ri 

and  its  variation  with  (rn/r^)  is  also  given  in  Table  V. 

piov  say  a  tunnel  with  a  radius  of  ten  feet  were  used  as  a  shelter. 
Assume,  toe,  that  the  boundary  is  soil-backed  and  therefore  the  maximum 
outermost  nodal  point  l*  ten  feet  past  the  outer  surface  of  the  boundary. 


1  v  4  Ti  -  Tn 
q  -  *  k  a - 

0  1  iir 


heat  conductance  area  at  the  average  radius 
\2 

ir.  +  r 


(r,  ♦  O 


5  ♦  {  £>  ♦  i  (5s) 


(4) 

(5) 

(6) 
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Figure  E-l  Nodel  Point  Rsdii  in  •  Cylindrical  Shelter 
Used  in  Determining  the Error 


1 

1 

] 

i 


«  VI  * 


TABLE  V 


♦- 


ERRORS  OF  C0NSIE2RING  THREE-DIMENSIONAL 
HEAT  TP ANSWER  AS  ONE- DIMENSIONAL 


r 


n 

II 

0(Cyli  ndrical) 

0 (Spherical) 

1.0 

1.0 

1.0 

1.1 

1.001 

1.002 

1.2 

1.003 

1.008 

1-3 

1.006 

1.017 

1.4 

1.010 

1.029 

1.5 

1.014 

1.042 

1.6 

1.018 

1.056 

1.7 

1.023 

I.O72 

1.8 

1.029 

1.089 

1.9 

1.034 

1.107 

2.0 

1.040 

1.125 

2.1 

1.046 

1.144 

2.2 

1.051 

1.164 

2-3 

1.057 

1.184 

2.4 

I.O63 

1.204 

2-5 

1.068 

1.225 

niCAN  nuirAficn  wtr-iON 
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Therefore, 


r1  +AP  + 


10 


(e: 


where  p  is  the  thickness  of  the  boundary.  A  reasonable  value  of  p 
would  be  3  feet  and  thus  the  error  in  the  heat  flow  due  to  u*ing  a  one- 
dimensional  analysis  method  of  a  three-uimensional  analysis  is  using 
the  cylindrical  surface  analysis  or  l8.4£  using  the  spherical  surface 
analysis.  (See  Table  V.)  Generally,  a  tunnel  designed  as  a  mass  shelter 
can  be  expected  to  have  a  radius  greater  than  ten  feet;  therefore,  the 
error  introduced  would  be  correspondingly  less.  Hence  a  shelter  boundary 
that  is  curvilinear  can  be  treated  as  a  flat  surface  with  an  error  given 
by  equation  (6)  or  (7). 
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ATPENDLC  F 

ESTIMATION  OF  ERROi-,  ASSOCIATED  WITH  CALCUIATION 
_ OF  AIK-CONDITIONING  LOAD _ 

The  performance  of  an  air-conditioning  system  is  a  function  of  the 
environmental  conditions,  thu  psychrometric  ctate  of  the  air  entering 
the  unit,  the  performance  characteristics  of  the  components  comprising 
the  system,  the  flow  field  in  the  unit  and  the  characteristics  of  the 
heat  sink  and  source  reservoirs.  Procedures  are  required  to  define  an 
air-conditioning  system  for  primarily  cooling  the  air  introduced  into 
the  shelter.  The  operating  principles  for  such  a  unit  can  be  based  upon 
a  large  number  of  choices.  However,  in  almost  every  choice,  a  coil  will 
be  utilised  for  cooling  and  dehumidlfying  the  air.  Theiefore,  the  SHEP 
code  has  been  constructed  so  as  to  require  only  the  specification  of  the 
coil.  It  is  assumed  that  once  the  coil  has  been  defined,  the  designer 
can  develop  a  system  to  produce  the  specified  coil  performance. 

In  particular,  the  coil  is  specified  by  denoting 

1)  a  coil  bypass  factor, 

2)  an  effective  coil  surface  temperature,  ind 

3)  an  air  flow  rate  across  the  coil. 

It  is  known  that  the  coil  bypass  factor  and  effective  coil  surface  tem¬ 
perature  are  a  function  of  the  thermal  load,  i.e.,  the  psychrometric 
state  and  flow  rate  of  the  air  through  the  coil.  In  the  SHEP  code,  this 
load  dependence  is  neglected  (since  performance  is  a  function  of  the  design 
of  the  entire  system)  and  the  bypass  factor  And  effective  coil  surface 
temperature  arc  assumed  constant.* 

•This  procedure  is  a  generalisation  of  a  coil  selection  procedure  presented 
in  "Fundamentals  of  Psychror-ctri cs  -  Fnrt  It",  1200-20,  Carrier  Air  Con¬ 
dition?  ng  Company,  page  19. 


OCNCRAi  AMERICAN  REP6ARCH  DIVISION 

F-l 


The  accuracy  of  this  procedure  has  been  checked  by  comparison 
with  the  performance  data  of  an  actual  air-conditioning  unit,  see 
Table  VT  ,  It  was  found  that  neither  the  deduced  bypass  factor  nor 
the  deduced  effective  coil  surface  temperature  are  constant.  However, 
if  median  values  of  the  range  of  each  of  these  parameters  are  used  to 
calculate  the  capacity  of  the  unit,  this  calculated  unit  capacity 
differs  from  the  actual  unit  capacity  by  at  most  37$  and  generally 
by  less  than  12$.  This  agreement  is  considered  sufficient  fcr  the 
inclusion  of  this  simplified  specification  for  the  air-conditioning 
equipment  in  the  SHEP  code. 
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APPENDIX  G 

TYPICAL  THERMAL  PROPERTIES  AND  DENSITIES  OF  BUILDING 
MATERIALS,  SOILS,  ETC. 

BUILDING  MATERIAL 


Thermal 

Conductivity 

Specific 

Heat 

Density 

CON 

(Btu/hr-f t-°F) 

CP 

(Btu/lb-°F) 

ROW  . 
(lb/ft3) 

Building  Board 

1. 

Asbestos -cement  Board 

0.33 

0.270 

120.00 

2. 

Ml 

Gypsum  or  Plastic  Board  4 

C.10 

0.270 

50.00 

M  „  „  1" 

2 

0.09 

0.270 

50.00 

3- 

Plywood 

oTo7 

0.5C0 

34.00 

4. 

Sheathing,  Wood  Fiber 

0.03 

0.500 

20.00 

5. 

Wood  Fiber  Board 

0.03 

0.500 

26.00 

6. 

Wood  Fiber,  hardboard  type 

0.12 

0.500 

65-00 

Finish  Flooring  Materials 

1. 

1" 

Cork  tile  -  g 

0.03 

0.435 

15.60 

2. 

Terrazzo  -  1" 

1.04 

0.250 

150.00 

Insulating  Materials 

A. 

B’-nkct  &  Batt 

1.  Cotton  Fiber 

0.02 

0.362 

0.8 -2.0 

2.  Mineral  wool,  fibrous 
form  processed  from 
rock,  slag,  or  glass 

0.03 

0.30 

0.5 

0.02 

0.30 

1. 5*4.0 

3-  Wood  Pi  be** 

0.02 

0.30 

3- 2-3-6 

CtNtMl  tMCmr^N  BfPtARCK  DIVISION 

G-I 


4- 


IV 


CON 

(Btu/hr^ft-6F) 

CP 

(Btu7lb-*F) 

ROW  , 
(lE7?t3) 

B. 

Board  &  Slabs 

1. 

Cellular  Class 

0.03 

0.185 

3 

2. 

Corkboard 

0.02 

0.485 

6. 5-8.0 

0.02 

12  _ 

3- 

Glass  Fiber 

0.02 

0.185 

4-9 

4. 

Expanded  rubber  (rigid) 

0.02 

0.339 

4.5 

5- 

Expanded  polyurethane 
(thickness  1  in.  and 
greater) 

0.02 

0.20-0.25 

1.5-2. 5 

6. 

Expanded  polystyrene 
extruded 

0.02 

0.30-0.35 

1-9 

molded  beads 

0.02 

0.30-0.35 

•  1.0 

7. 

Mineral  Wood  with 
resin  binder 

0.02 

O.30 

15 

C. 

Loose  Fill 

1. 

Mineral  wool 
(glass,  slag  or  rock) 

0.03 

0.30 

2. 0-5.0 

2. 

Perlite  (expanded) 

0.03 

0.378 

5. 0-8.0 

3- 

Sawdust  or  shavings 

0.04 

0.30 

0.8-15 

4. 

Silica  Aerogel 

0.01 

0.205 

7-6 

5. 

Vermiculite  (expanded) 

0.04 

0.04 

0  0 
•  • 

U*  U* 

O  O 

7. 0-8. 2 
4. 0-6.0 

6. 

Wood  Fiber,  redwood, 
hemlock,  or  fir 

0.02 

0.30 

2. 0-3. 5 

Masonry 

Material* 

1. 

Cement  Mortar 

0.42 

0.271 

116.  C 

2. 

Gypaur.- fiber  concrete 

87-1/2,*  gypsum, 
l?-l/«**  wood  chip* 

0.0* 

0.25 

51.0 

QCMFtlAl  AMERICAN  «Cr*»  M)'H  ©IV»y*C?i 

G-? 


•  1 
-  J 

.'3 


:: 

CON 

(Btu/hr-f t-°F) 

CP 

(Btu7Tb-°p) 

ROW  , 

(i57n3) 

3. 

Lightweight  aggregates 

0.43 

0.20-0.25 

1 

including:  expanded  shale, 

120 

J 

clay,  or  slate;  expanded 
slags;  cinders;  pumice; 

0.30 

0.20-0.25 

100 

3 

perlite;  veraiculite;  also 
cellular  concretes 

0.21 

0.20-0.25 

80 

i 

0.14 

0.20-0.25 

60 

J 

0.10 

0.20-0.25 

40 

3 

0.08 

0.20-0.25 

30 

0.06 

0.20-0.25 

20 

3 

Stucco 

0.42 

0.20-0.25 

116 

V  Masonry  Units 

1. 

Brick,  common 

0.42  ‘ 

0.20' 

120 

3 

Brick,  face 

0.75 

0.20 

130 

3- 

*1 

Clay  tile,  hollow 

J 

1  cell  deep  3  in. 

0.31 

0.15 

30 

_ 1 

2  cells  deep  6  in. 

0-33 

0.15 

30 

■n 

2  cells  deep  8  in. 

0.36 

0.15 

30 

3 

2  cells  deep  10  in. 

0.38 

0.15 

30 

i 

3  cells  deep  12  in. 

0.4o 

0.15 

30 

j  !•. 

] 

Concrete  blocks, 
rectangular  core 

Sand  and  gravel  aggregate 

0.15 

0.25 

30 

J  5. 

Stone,  lime  or  g*nd 

1.04 

0.22 

150 

J  VI  Plas 

Urir5  Materials 

1. 

Cerent  Plaster, 

J 

sand  aggregate 

0.4? 

0.20 

H6 

«.  . 

Cyprus  Plaster 

% 

J 

} 

Perlite  aggregate 

Sand  aggregate 

Versucui'*e  assregaie 

CM? 

0.47 

0.01 

odd 

45 

105 

4> 

J 

DENEr.Ai  AI.'iOCAN 

esr-EArtCM  Oivtc 

0-3 


VII  Roofing 


CON  CP  ROW  - 

(Btu/Hr^ft-'P)  (Btu7Tb-*F)  (li>7Ft3) 


1. 

Asbestos-cement  shingles  0.40 

0.20 

120 

2. 

Asphalt  roll  roofing  0.5U 

0.20 

70 

3. 

Asphalt  shingles  0.19 

0.20 

70 

VIII 

Siding  Materials  (on  flat  surfaces) 

1. 

Shingles  -  asbestos-cement  0.40 

0.20 

120 

2. 

Siding 

1" 

asbestos-cement  -  ^  0.10 

0.20 

120 

asphalt  roll  siding  0-.54 

0.20 

120 

3- 

Architectural  glass  0.83 

0.15-0.20 

150-175 

IX 

Woods 

1. 

Maple,  oak,  and  similar 

hardwoods  0.09 

0.50 

45 

2. 

Fir,  pine  and  similar 

softwoods  0.07 

0.50 

32 

CONFIGURATION 

Thermal  Conductivity 
Unit  Distance 

Specific 

Heat 

Densi  k.y 

,  m/M 

(Btu/hr-f t-*F) 

CP 

(Btu/lb-T) 

ROW  , 
(lb/ft3) 

I 

Vertical  Space  1.03 

0.24 

0.075 

II 

Kori rental  Space 

1. 

Heat  flow  up  1.18 

0.24 

0.075 

2. 

Heat  flow  down  1.01 

0.?4 

O.C75 

lit 

!»5* 

In  Mncd  Space 

1. 

Heat  flow  up  1.11 

0.24 

0.0^5 

2. 

Heat  ft  -'vdevr.  1*13 

0.24 

0.075 

AWrniCAN  Brf'EArCH  OtVIsEIOM 


SOIL* 


Thermal 

Conductivity 

Specific 

Heat 

Density 

con 

(Bt"./hr-ft-°F) 

CP 

(Btu7lb-4 

F) 

ROW  , 

(W‘t3) 

I 

Crushed  quartz 

0.96 

0.188 

100.00 

II 

Graded  sand 

O.83 

0.210 

100 -CO 

III 

Gravel 

0.75 

0.200 

no.  oc 

IV 

Crushed  feldspar 

0.50 

0.190 

100.00 

V 

Crushed  granite 

0.46 

0,230 

100.00 

VI 

Sandy  loam 

0.54 

0.210 

110.00 

vi:.; 

Crushed  trap  rock 

0.42 

0.195 

100.00 

VIII 

Sand 

0.37 

0.210 

100.00 

IX 

Clay 

0.33 

0.220 

100.00 

X 

Silt  ioan 

0.42 

0.220 

90.00 

XI 

Silty  clay  loam 

0.42 

0.200 

90.00 

BOUNDARY 

ORIENTATION 

Direction  of 
Heat  Flow 

Film  Heat-Transfer 
Coefficient 

HI,  HQ 
(Btu/hr-ft 

-T) 

I 

Vertical 

Into  or  out  1 
shelter 

■jt 

1.46 

II 

Her i run ‘ t i  (ceiling) 

Into  shelter 

1.08 

Out  of  shelter 

1.63 

III 

Kof  I  ( fie  >r } 

Into  shelter 

1.63 

Out  of  shclt< 

?r 

1.08 

Ofurn.M.  dMM»;c»N  r,f.T:«rcM  oivir.in%‘ 


3 

3 

3 

') 

j 


j 


i 

J 
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APPENDIX  H 

SHEP  CODE,  MODIFICATION  3 
FLOW  CHART  ANI)  LISTING 


O'.'Ni 


tf  l  r.r.- 


:a\ 


nr»n 


4 


PROGRAM  SHEP 
C 

C  SHELTER  ENVIRONMENTAL  PREDICTION  CODE  ISHE®— MOD.  3) 

C  GENERAL  AMERICAN  RESEARCH  DIVISION 
C  GENERAL  AMERICAN  TRANSPO® 7a f ION  CORPORATION 
C 

C  THIS  COOE  CALCULATES  THE  RESPONSE  OP  A  SHELTER  TO  TIME-VARYING  THERMAL  LOADS 
C  AND  VENTILATION  INLEI  PSYCHROMETR1C  CONDITIONS.  THE  CODE  ACCOMMODATES  SOLAR 
C  RADIATION.  BOUNDARY  SURFACE  HEAT  TRANSFER.  MOISTURE  CONDENSATION  ON  THE 
C  BOUNDARIES'  EQUIPMENT  AND  LIGHTING  LOADS*  AND  AIR-CONDITIONING.  THE  USER 
C  MANUAL  EXPLAINS  THE  VARIOUS  FEATURES  OF  THE  CODE.  PRESENTS  THE  INPUT  AND 
C  OUTPUT  FORMATS.  GIVES  REPRESENTATIVE  INPUT  DATA  VALUES.  AND  INCLUDES  AN 
C  APPLICATION  OF  THE  CODE  TO  A  TYPICAL  SHELTER  STRUCTURE. 

C 

COMMON/MATRIX /A  I NV  t  40.40.20  I 

DIMENSION  TWALLI 40 *20)»A(40»40)«I NDEX 140.2 ) »T 1 40) 

DIMENSION  EOT IMl 365 1.T02PI 336 ).TWBI 336»  »TEF? 13361 

DIMENSION  PSATU20)  .INI  5. 20  (.DPI  5.201. DX<  5 .20)  .ALPHA!  5.201. CONI  5. 2 
10) .ROWI5.20) .CPI5.20I »CW<20) 

DIMENSION  INT ( 60 ) .OSi 2 .20) «BB( 40 ) »BMUL 1 40 *20 ) .TDI I  30 ) »TWl I  30 ) • TEI I 
ISO)  .T0IAV13C)  .TWl AVI 30)  .TEI AVI  30)  .TSOI LI 20) 

DIMENSION  XI (20) .A2I20) .LAVMAX 1 20 > .OLATCI 20 > .HI  1 20 ) .HOI 20) .NMAXI20 
1 ) . I TMALL 1 20) .0SWI20) .owl  20) .AREA  1 20 ) . TDDI 14 ) .T wDl 14) • TED  1 14 ) • TDAvG 
2  1 14)  .TWAVGI 14 )  .TEAVCil  14  )  •  10ATA.  1 13  I  »L5T  1 2 )  .TRANS « 20  ) .  ASOILl  20 ) . 

3AJITI20) 

JJJJ-1 

INITIALIZE  ALL  MATRIX  ELEMENTS  TO  ZERO. 

1311  00  1  1-1.40 
DO  1  J-1.40 
DO  1  X*1 .20 

1  AINVC I.J.KI-0.0 
DO  2  1-1.40 
TID-O.O 

DO  2  J-1.2 

2  INDEX  1 1 .  J) *0 
DO  X  1-1.40 
00  3  J- 1 *20 
BMULl I.JI-0.0 

3  TWALLI I . JI-0.0 
00  5  l-l .365 

3  COT IMI) 1-0.0 
DO  6  1-1.336 

4  TD2Pm-T*0l  I  l-TEFFII  J-0.0 
DO  7  1-1.120 

7  PSATI  D-0.0 
DO  6  1-1.5 
DO  •  J- 1 .20 
IN  1 1 »J)-0 

«  OPII . JI-OXI l . Jl-ALPMAI I .J 1-CON  I  I .J) -ROW  1 1 *jl-CPI I.JI-0.0 
DO  «  1-1.60 
t  inti:— o 

DO  10  1-1.2 
DO  10  J-1.20 
10  QS<I.J>*0.0 
DO  302  1-1.30 

302  TDI II I-TWII ) )*TEH I l-TOIAVl I l-TwlAVf I I-TEIAVI! )-0.0 
DO  30)  1*1.20 
vWI I  1*0.0 


♦ 


LAYMAXI  |  )*NMAX  <  1 l-ITWALLI I  )*0 

)03  XI U  >-A2(l l-OLATCt 1 >*H! 1 1 >*’lOl I  )*Q$W( I ) *Q*M I ) ■ AREA ( I )*TRAN$( l )*ASO 

-  llLUJ*tSQlLll)»AJIILIl*0.0  . 

00  3U4  1*1*14 

304  TDOtll-TMIlWEDI  i>TDAVG<  1  >*TWAVGt  I  )*TEAVGI  1 1-0.0 
00  303  1*1*13 
303  {DATA! 1 )  *Q 
00  306  1*1*2 

306  LStm-0  . . . . 

C 

C  READ  INPUT  DATA.  AND  SET  CONSTANTS* 

C 

IP (JJJJ-1 1 307*307*308 

307  00  11  1*1*363 

U  REA0-1001»£QTIM1U . .  .  .  .  .  . 

00  12  1*30*120 

12  READ  1001  * PS ATI! I  . 

308  READ  1006* IDATA 

REAO  1003«ALT»ALAT*ALONG»NW*GC'*HH»W»P»F1*OUITE*OEQUIP*0»MRS*LST(1) 
l»NO»NT 

_ 00-14  J*  1 «  mw _ 

READ  l020*AREAI J)  *LAYMAX( J) .NMAXI J) *Ml < J> *H0( J) *AS0IL(J) 

IF (ASQILtJ)J 1312*13*13 
1312  READ  1001 «TR ANSI J) (XI I J ) 

IF  (XlUL-90.0)  14*13*999 
14  REAO  1003  *AZU) 

13  L£M*JLAYMAXU  1 . . 

00  14  K*1 *LEM 

16  HEAD  -1009»DPlX»J)»OXtX* J) «C0N(X*J) »ROW(X* J) *CP(K»J) 

READ  1013*CFHIN*TC01U»WC01C*6F*F2  *T02*R2  »TWB ( 1 ) iTEFF 1 1 ) 

_ GO  18.J-1.NW . . 

1FIASOILIU) >18*18*1703 

1703  -READ  -lCQSxTSOlLU) _  .  .  .  .  . 

18  CONTINUE 

IF UQATAI1 ) 119*19*21 

19  READ  1008 *TO 

DO  20.J*1»NW  . . 

NEM*NMAX ( J ) 61 

- -DOaQ_JL«l*iiEJl _ _  _ _.  .  . 

20  TWALL(K*U)*TO 
GO  TO  23 

21  DO  22  U*1*NW 

NEK*NMAXIJ)61..  ....  . . . 

DO  22  K*1*NEM 

—  ftEAU.-10.QJ*  t£L _  . 

22  TMALLUiU)*TO 

23  REAO  1001*0.1  . 

’MR5*MR$ 
rtODAYS*!MRS/24 
1FINOINTI 1701 *1701* 1702 

1702.  JlQ0»2?UQlia . .  .  . 

00  17  I  •  l »  NOO 

17  REAO  1006»INTI1) 

1701  10*0 

LHN*2S 
R1N0*0 

OLM*>0«.0 _ 

AOT-O.O 
tNRNO»i 
NTAU“0 


PI -3. 1*159265 
OM-PI/12.0 

*■1.0 
MONO-1 
6IL-TD2 
V5»6C*HM»W 
READ  1021 »TA1 .TA2.T01 1.TD12.R11.R12 

CALCULATE  INTENSITY  OF  SOLA*  RADIATION  FOR  EACH  EXPOSED  BOUNDARY . 

25  LST ( 2 J -LST ( 1 IS1 
DO  53  I-1.2 
LSTM«L$T(II 

251  DO  53  Jal #NW 
IFUSOILUM252.53.53 

252  IF ( IDATAC 10) 126*25.27 
'  26  C2-294. 

C3-3.183 
60  TO  28 
2?  C2-200. 

C3-2.228 

28  IFUIUI-90.)  29.29.999 

29  ETAR-0.<»10152*SINF  1 1 1 0*80 )*P  1/182* 5  ) 

S1«TANF(£TAR) 

52- SINFIETARI 

53- C05F ( ETAR ) 

PH1R-ALAT/0N 
S6-TANFIPH1RI 
S5-SINF (PhIR ) 

S6*C0SF i PHI* ) 

A2R-AZUI/0N 
XIR-XI  UI/ON 
S7»C0SF«XIRI 

58- SINF  tXlR) 

IFIPHl*-(Pl*O.Sn  30.998.998 

30  IMETARfcPHlR-IPI*O.S)l  31.998.998 

31  IFIETARI  34.34.32 

32  IFtPMlR-IPI*0.5-ETAR})  34.33.33 

33  ThETAR-0.0 
60  TO  35 

34  ThETA*-AC0SF(S4«SM 

33  ThETaS-2.*P!-TH|TAR 
IN* ALONG/ 15 

Tl-LSTN4E0T!M( IDl-ALE  5.041N 
TmETA»0m*TI 

59- C0SF i THETA ) 

IFITMETA-THETARI  36.37.37 

34  DIR-0.0 
01F-0.0 
60  TO  52 

37  JMThETA-ThCTAS!  38.38.36 
36  BE  TAR -AS INF J SS*S2~S4*S9-S3> 

SlO-COSF IBETAR I 
SU-SINF  I  BE  TAR  I 
IFIBETARI  36*39.39 
39  |FteCTAR-tPi*0'»M  392.391.392 
391  6AK.VAR-AAR 
60  TO  393 


392  GAMMAR«AC0SFJl»/S10*«S6*S26S5*S9*S3n 

393  IF(GAMHAR)  40.40.41 
40  GAMMAR*GAMMAR6P  t 

4J  DIRNORN*C2*U»-EXP‘-C3*B£TARn 
S13"S10*C0SF 1GAMMAR-A2R 1 
AUPHAR>ACO&F  ( S13*S76S1 1*S6 ) 

$12*COSF t ALPHAR) 

IFULPHAR-<Pl*0.3>)  43.43*42 

42  DIR*0.0 
GO  TO  49 

43  IFIS13I  43.44.44 

44  FALPHA*0.35*5i3**2.60.45«Si360.5S 
GO  TO  40 

43  IF IS13&0.04)  47.47.46 

46  FALPHA*0. 625 *$13**2.60.523*5 1360 *53 
GO  TO  40 

47  FAWPMA*0. ’ 4 

40  DIR*DIRN0RM*S12 
49  OIF1-$11*(FAUPMA*(1.-S0)6S0I 

OIF2»0.677E“O6*DJRNORM**2-l*5E-O4*DIRNORM60.01#3 
01F*01F1/DIF2 
32  T0T*0IR6DIF 
QSU  .JI«TOT 
53  CONTINUE 
OST-O.O 
OSOUARaO.O 

DO  54  J*1«NW  - 

OSWl J)«OS(2. J) 

QSOL AR-TR ANS t  J I *OSW ( J 1 • AREA ( J ) 6O0OLAR 

34  O0T*OST6OSl2.J)*AREAt J) 

FOR  EACH  SOUR  OF  OCCUPANCY  AFTER  THE  FIRST.  RESET  HOUR  NUMBER  AND  SET 
CONDITIONS  AT  The  CNO  OF  The  PRECEDING  hour  TO  THOSE  AT  THE  BEGINNING 
OF  THE  SUCCEEDING  HOUR.  READ  THE  INUET  CONDITIONS  FOR  THE  SUCCEEDING  HOUR. 

IF ( IHRNO-1 157.24.37 

35  DO  36  J*1.NW 

56  OSU  .J»«QS(2.JI 
LSTM-CSTMtl 
IHPNO* 1 HRN06 1 
IF ( 24-LSTMI3612 *5612 .251 
3612  tSTM-0 
1D*ID61 

IF  1 365" 10 1 3611 .251*231 
3611  tO-1 

GO  TO  231 
37  TD1 1*T012 
Rll*R12 
TA1*TA| 

AOTaOT 

IF  1 1 DAT At  121 1111. 161. 112 

161  READ  1016.P.F1.CFMIN.OUTC*QCOu|P 

162  IFUOATAUJH16J.16J.60 

163  REAO  1019.T012.TA2.R12 
GO  TO  60 

36  03T*0.0 
OSOCARaO.O 
DO  59  ,-l.NU 

OSO*OSa.J>6ADT*tOS(2.J>«0«U*J>  » 

0S«< JI.OSO 


* 


OSOLAR* TRANS ( J ) *OSW ( J)*AREA( J) 4QSOLAR 
59  0ST-0ST60S0 
C 

C  FOR  TIME  INCREMENTS  DURING  THE  HOUR*  DETERMINE  INLET  CONDITIONS  AS  A 
C  LINEAR  INTERPOLATION  BETWEEN  HOURLY  VALUES* 

C 

*0  TD1»TD11-ADT*ITD11-TD12I 
R1*R11-A0T*CR11-R12» 

PG*5*132*EXP|0*0329*TDil 

PGPP*Rl*PG 

8P-0.9B£-04*ALT**2-0.07S9*ALT42114.2 

m*IO*422*PGPPI/IB?-PGPPI 

TA*TA1-ADT*ITA1-TA2> 

LL-1 

C 

C  CALCULATE  NEw  tNLET  CONDITIONS  WITH  AIR-CONDITIONING*  DETERMINE  PROPERTIES 

c  op  the  inlet  air  for  the  beginning  of  the  time  increment* 
c 

I F 1 1 DAT  A ( 7  > >403**03*402 
402  TD1-TD2MF2/F1IMTD1-TD2I 

PGPP«S.132*R2*EXP»0.0329*TD2I 
W2" I 0.622*PGPP 1 / (BP-PGPP i 
W1*W24IF2/F1 »*IWl-W2) 

ADMIX*  CF1-CFMIN»/FHBF»U.-!(F1-CFM!N! /FlU 
7DA»TC01L4ADMIX*ITD1-TC01L) 

WAaWCOI L4A0M I X* I Wl-WCOt  L ) 

Ei*EXPI0.0329*TDA)  --  ~  - 

PGPP>WA*BP/(WA40*422I 

PG*5.132*E1 

R1«PGPP/PG 

WACa40.0*CFMlN*0T*I0.0lB744*BP-0*0942*Al*EimTDA44S9.491 

HF«i093*9-0,»b7«TCOIL 

PAlRCON0a0.24*WAC*TTDl-TDAJ/0T4WAC*tWi-WAI*HF/DT 

TD1-TDA 

W1*WA 

GO  TO  SI 

401  paircond-o.o 
SI  ElaCxP(0. 0329*1011 
E2a£XP(0. 0329*1021 
HA$1«0«24*T02 
HVSla104l. 040. 444*102 
HA1«0*24*TD1 
KVla 1041. 040.444*101 
MFGal09S*9“0*547*TD2 

WASl*VS*(0.01l74b*BP-0*0942*R2*E2l/(T02S439.69l 

WVS1"VS*(0*039«7*R2«I2I/(T02S439.491 

WA1«40*0*F1»0T«(0*01S744*BP-0.09S2*R1*C1 1/11014439*491 

WVl«40*0»Fl*OT*tO.OS9l7«Rl*Ell/<T01S439*49l 

HS1«WAS1 »HAS14WVS1*HVS1 

H1*WA1«KA14WV1«HV1 

WAO*WAlwWASl 

C 

C  CALCULATE  Tnr  SCAT  TRANSFER  through  The  boundaries  out  TO  temperature 
C  DIFFERENCES  AND  OWE  TO  MOISTURE  CONDENSATION. 

C 

IF  t I0ATAI41 149*49*70 
49  OWC-O.O 
OLH-C.O 
NT AJ-NT4W41 

>H>N*lHR*iO~l  . 

H'7 


GO  TO  146 

70  00  71  J-i.NW 

ILL-0 _  _  . 

LEM-LAYMAX(J) 

00  71  K-l.LEM 

lN(K»J)aDP(K.J)/OX(K. J)G16ILL 

71  ILLMNiUtJ) 

JW  ALL-1 

WAY-1  .  _ 

I -l 

72  MCN*AMAX t  JWALU61 
00  4  tNAai.fl£M 
DO  4  JNA'l.AEM 

Al INA.JNA) a0.0 
4  BB  UNA) -0,0  _ 

72 1  ALPHA ( WAY JW ALL I -CON (LAY  , JWALL )  /  ( ROW t LAY  »  JWALL )  *CP  ( LAY  i JWALL I ) 

BAS! L* 1*06 CON I LA Yt 1  * JWALL l*DXt LA Y  * JWALL I / ( OX ( LA YG 1  .JWALL I *CON ( uAY * 
1 JWALL  ) 1 GROW ( LAY  •  JWALL  t  *CP ( LAY  »  JWALL I -DX ( LAY  «  JWALL ) **2/ 1 2 -O-CONl LAY 
2 • JWALL I *07  I  GROW '.LAY  61 * JWALL  >  *CP ( LAY  61 • JWA ' L ) -OX ( LAY » JWALL ) *DX ( LAY 6 
91 .JWALL) /(2.0-CON (LAY  » JwALL )*DT ) 

7  9 -XrilalUbiltyx&l _ 

74  !702>T02 

PS-PSAT I ITD2 ) G( TD2- 1 T  02 ) * ( P$AT ( I T02  61 > -P5AT ( IT02) ) 

P0-P6-R2 

1 T WALL ( JWALL ) -7 WALL ( 1 1 JWALL  ) 

1 CLP 1 * 1 T WALL ( JWALL  )- 

.  CL  F  2"  tWALL  UutJw  A  LL 1  _  . 

P4W-P5AT  ( ! ELF 1 ) 6 ( ELF2-IELF 1 ) • ( P5AT 1 1ELP161 )  -PSA?  1 1  ELF  1 )  I 
TENV" IID26TWALLt 1* JWALL) 1/2*0 
r.OA)R-9.29E-07-TENV— 2.0-1.46E-04-TEWV60.0969 
7 7  C»2 :i 9249-HFG-Hl  ( JWALL J • f  PO/ ( BP-PO )-P5W/  ( BP-PSW >  I  /RQA l R 
CWt JWALL) -C*AREA(JWALLI«OT/HPG 

- lPLLCwlJWALUJ.2I2t79.JU  ...  .. 

771  AJ1 T ( JWALL I-AJIT (JWALL I 6CW( JWALL) 

GO. TO. 79.. 

771  IF (A65P (CW( JWALL  * l-AJIT « JUALL)) 771.771.774 
774  CW(JWALL)— AJ1TIJWALU 

C-CW( JWALL >  *HFG/ (AREA  t JWALL ) -DT ) 

-GCLTO-IU _ _  _ _ 

79  OLATC( JWALL) >C 
IF (NTAUl9O.Q0.tl 

•0  A(I»n-i*OfrCON(LAY«JWALL)/(Hl(  JWALL  )*DX(  LAY  .JWALL))  GROW  (LAY*  JWALL) 
l*CP ( LAY .JWALL ) *0X( LAY .JWALL ) / ( 2.0-H I ( JWALL ) *0T ) 

A( I » Itl ) «-CON ( LAY. JWALL )/(Hl ( JWALL) *0X1 LAY t JWALL ) ) 

•1  6RU I  -ROW  l  WAY » ALL )  •  CP  ( L AY  t  JW ALL  I  »0X  ( LAY .  JW ALL )  /  ( 2. 0»H  1 1  JWALL )  *07 
1 )*TWALL( 1 .JWALL ) 6T02GL/HJ (JWALL) 

1-161 

it  IP ( 1—1N( LAY. JWALL) )99.t9.t6 
19  IP (NT AU) 84 .94.99 
14  All.t-1)*- 1.0 

A(  1 1 U POKCLAlY  »  JPALL.)p*2/(  ALPHA  (LAY.  JWALL )  *07  )62  * 

Atl.IGl)— 1.0 

99  BB ( 1 ) *0X ( LAY .JWALL ) -2- T WALL ( I .JWALL ) / ( ALPHA ( LAY .JWALL ) -07 ) 

I-1G1 
GO  TO  92 

•6  IP (LAY**LAYHAX(  JWALL)  197.90.996 
97  IF (NTAUiHtOfl i9.«  _.  _ 

99  All. I -1)— 1.0 
A ( | *  t I-BA5IL 

A ( I . I G1 ) —CON ( LAY4 1 . JVALL ) -OX ( LAY .JWALL I  /  ( OX ( L AVO 1 »  JWALL ) -CON I LAY  * 
1 JUALL)) 

BO  C3( I )-( ROW ( LAY  *  JWALL > -CP ( L AY  * JWALL ) »0X ( LAY . JWALL ) « -2  *6R0W ( LAY 6 1 . JW 

Hi 


S; - - - 

2/12*0 *CON (  LAY » JWALL » *0?  > 

IF t LAY-LAY MAX  (  JVALL 

,0  SliSlLUHAUllWWW 

91  TAJ-TA 
60  TO  93 

92  fAJ*0*5#<T02ATAI 

"  "iH3S:;:sss:si!s:siss::» - - 

•s' ir  usoiil*  1  iwi  •*»>•** 

. . . 

jALLl I*TwALL( I • JWALL  > 

??,ISJiJ«jwALL»-moe»ioo»m 

. ~ . . . . . 

i»*on  . . 

|MNTAU»101*10-  *102  - 

101  A(t»l-l,“i*°  ,.y 

103  lF(NTAU»104tl04«13Si 
10i*  MtM*NMAX« JWALUH 
00  10t  !•! »NEM 
JOI  lHO€Xlt»l»"0 
!1«0 

109  AMAX— l* 

OA  110  1-ltNCH 
IF ( tHOCX ( |*1)>1 10 *111*1*® 
tit  112  J*l 

|FUNO€*<  J»i  niiJ»il5,Til 

113  T£MP*ABSr<Al  I  »->>  * 

|  F 1 1 CM?" AMAX > 1 1 2 • 1 1 2  * H 

H%  | ROW*  I 
ICO'-  •  J 
amax*»temr 
111  CO*TiWl 

»w  srjas„..»»».n» 

133  FOAHATtllH  2CR0  RIV0T1 

CO  TO  A* 

114  tKOE*UCOL»U»t*<»* 

|F  1 1 ROW- 1  COL  1 1 ^9 »llt  *11* 

119  00  UO  J"l»*£* 
t|MR*M  *ROw»>it 

At  JROW* J>*At ICOL»J' 

120  AUCOL*J»«TCw 

u*iim 

l«Ot** 1 1 

lit  RIVOT**UCOL»:  *.0L*  J #  ^ 

aucol*kol>-i»o  at  rn 

R I  VO  *  *1»/^1  VOT  f 


A  AAA  AAA 


*• 


DO  121  J>1*N£H 

121  A I  !COL*J)*A( ICOUfJ MP IVOT 
DO  122  t*l*NEM 
iFU-lCOL>123. 122*123 

12.3  TEMf»"AII«ICOL> 

All  •  I  COD  *0*0 
DO  124  J-1*NEM 

124  AI1»J)*A(» 1 J )»A( ICOL * V) »TEMP 

122  CONTINUE 
60  TO  109 

129  1 COL" INDEX  1 11*2) 

1ROM-INOEXIICOL*1) 

DO  124  I*1»NCM 
TEMP*A( I.IftOHl 
All » IR0W)"AI 1  * 1C0LI 
124  All* ICOL 1 "TEMP 
11*11-1 

229  1FII I >129*127*129 
127  DO  139  1*1 *NEM 

BMUL 1 1 »  Jw  ALL } *BB 1 1 1 
DO  134  J*1 *NEM 
194  A1NV 1 1 *U* JWALL i »A 1 1  * J I 
1991  IFI JWALL-Nwl 143*142*999 
1-3  JWA!.L»JWALL01 
1*1 

LAY-1  --  *  - 

GO  TO  72 

142  DO  144  K'1*NW 

1*4  QtHK  )*Hl IX l*AREAlX)*l TD2-TWALL 11*41  I 
OWC-O.O 
OLH-O.O 
00  144  P-i.NW 
OwC-OWCAQWl 21 

144  QLM"QLm9GLATCU»*AREAIJI 

CALCULATE  METABOLIC  HEAT  ENERGIES* 

NTAU*NTAUtl 

tHRN*tHRNO-l 

1*4  IF  ( T|FF  UMRNO  1-9  7.0  >149*  149. 147 
147  IF  » TEFF  1 1 HRNOI-102*0 1 148  *149  *994 
14|  OT«-l,909«TEFFl  J  -t»NO)**2L249.7*TEFF  1  1MRN0I-10794.2 
GO  TO  191 

149  IF  I  Tiff (IMRNO >-40. 0)99*. 140. 140 
140  CT»— l»492"TCFFt 1 HANOI 49 14.0 
19}  JFIT02-90. 01992*142.192 
19)  Qtt*S«-O.Q4t79*T02«mi.929*TD2t523.0 
OSCN«OSENS*F 
«LA*OT-OSEN$ 

QL*OLA*A 

CALCULATE  THE  SMCLTEfl  ORY-OUL*  TEMPERATURE  AND  RELATIVE  HUMIDITY 
AT  t*E  t*0  Of  The  T|H€  INCREMENT. 

19)  IFIR-) .011 94 *199*199 
134  OL*OL*R 

44 1  NT  200 1.| man. n 

133  00»MiU.Ml4l05lN<-ttL*IMVSl/MFG)-OvC40L:iES*'7t0U|R-PAIRC0N04OS0LAR-0LH 
1  1*07 

MVL*0L*07  AmFG 

fcto 


1#VC0ND-QLH»DT /MFG 
WVO-HVIGWI  *  1-WVCOND 
IF IWVCOND) <  \554.1556 

1556  IF lwv0l 155>»  '"*1534 

ISM  WVO-v.O 

1554  TD2PllHRN01»*DQ-1061.0»WVO>/«0.444#WVC&0.24*WAO> 
IFIABSFIGIL-TD2PUMRN0)  1-1.00  i  1552  •  1551 .1551 
1552  LI* I DATA? 3 > 

GO  TO  163 

1551  IF (T02PI 1HRN0 I -50.0 1986.156. 156 

156  IF (WVO) 1561 .156 1.1563 

1561  PV2-0.0 
R2P-0.0 
GO  TO  1562 

1563  PV2*6P/I0.622»WAO/WV061.0) 

|P$-EXPFI0.0329»TD2P( I HANOI  I 
*2P*PV2/l5.132«fPS) 

1562  IFIR-0.01 157 .157 .158 

157  R2P-1.0 
GO  TO  K3 

151  1FIHONO-2) 159.160.160 

159  MONO-2 

IP (R2P-1. 01 163.163. 160 

160  IFIR2P-I.0il62.162.161 

161  A-R-DR 

GO  TO  153 

162  R-R60R 

DR-O.i-DR  ' 

MONO-MONO&l 
IF (MONO-4 1 161.161*163 

163  1FUL-IDATA<31 1164*165.991 

164  TD2*0.5*(T026T02PllHRNOn 
GIL-T02 

R2»0.5*(R26R2PI 

C 

C  CAwCULATt  ThC  SHELTER  HET-8UL6  ANO  EFFECTIVE  TEMPERATURES. 

C 

DO  1641  l»l.NW 
1641  AJlTItl-AJITIII-Cwm 
LL-LL61 
GO  TO  61 

165  DO  1555  X-i.N* 

NCM-KMAXtKItl 
00  1555  1-1. NCM 
TUI-0. 0 

OO  1555  J-l.NCM 
A(I.JI-AINV( I.J.XI 

bbui-bmuuj.ki 

Tm-TUI6AU.Ji*DD«J» 

1355  TVALUl.Kl-Tt’: 

IF (R2P-1.0I 167*166. 990 

166  TtfRt!HR*OI-TD2PUM*NOI 
tcffuhrno»*to»pi  imrnoi 
GO  TO  168 

167  T*o-T02Pi Imrnoi 
MC-10.0 

1 

168  PSS-5 .1 32 -CXPF 10.03 29*TVOI 

TW8( |HRNOI-( t  6*  -PSSI*T02PI lHRNO»-2«0O.O*IPSS-PV2> 1/1  #P 
1V2-2.3-OSSI 

tf'li 


6lO*P 


169  IFITWBI  lMRNOI-TWOU70.il  r?»171 

170  TWO«TWO-OWO 
GO  TO  169 

1 7 1  TWO'TWO&OWO 
DWO*0«1*DWO 
MM-MM61 

IFIMM-3U70.170.172 

172  TEFi«107.5*iT02PlIHRN0)-TWBI  !HRN0n662.3*TwBllMRN0) 

f  TCFF(lHRN0»»TEFi/?62.36T02P< IhRnO>-TwB( 1HRNO) t 

C  PRINT  RESULTS  FOR  THE  HOgR* 

C 

IPUIND-IDATAI5M 173. 174*969 

173  KIND-KINDS! 

GO  TO  180 

17*  IFU0ATAmni79.178.175 

175  IF ( JmRNO-I >988  *1 76* 177 

176  PRINT  2002 
PRINT  2037 
PRINT  2003 
PRINT  200* 

PRINT  2005 
PRINT  2030 
PR, NT  203i 

177  IH lHSNO-LMNI1772*r771*1772 

1771  LMN»LMN&2* 

PRINT  2038 

1772  F3^F1/P 

print  200*. ihrno.fs »td2pi  ihrno) , f ws < ihrnoi *teffumrno> .purcond.qw 

1C.0LH  iOSOLAR  tQSENtGL.  ,QL  I  TE  tCEQU'P 

178  IFI1HRN0-IHRSU79.187.987 

179  k:no«i 

183  T02»TU2P(l;iRNOI 
R2-R2P 
AOT-AOTSOT 

ItFFI IHRNC61 )®TEFF« IHRNO* 

2*  IF  <407-1*00001 158*55*35 

» 

» 

■  ERR  R  MESSAGES 

*99  PRINT  2007 * J.X I ( J) 

GO  TO  5000 
9*8  PRINT  200 i 
GO  TO  5000 
*97  PRINT  iviO* 

GO  TO  5000 
9’*  PRINT  2010 
r  *0  5000 
**5  T  2011 
OO  to  5000 

»*<  PRINT  2012*  |HR;  « T£FF ( IHRNOI 
GO  TO  1*5 

9*3  PRINT  2'12, I«RN*rEFFI IHRNOt 


GO  TO  1*5 
9*3  PRINT  2 012 
GO  TO  :  %■' 
9*2  PRINT  201* 
GO  TO  IT. 
»*1  PRINT  2Ci5 
GO  TO  1*5 
*90  PRINT  201* 
GO  TO  16* 


2  il*«  IHR(.»’02P»  IKTNOI 

\  *.  * 


nAn 


989  PRINT  2017  '%• 

60  TO  174 

986  PRINT  2018 

60  TO  176 

987  PRINT  2019 
60  TO  187 

986  PRINT  2036.TD2PI 1HRN0) »iHRN 
60  TO  196 


CALCULATE  ANO  PRINT  TIME-AVERAGE  SMELTER  TEMPERATURES. 


187  T0MAX-T02PI 1 J 
IHH" JHRS-1 
LM*1 

168  00  190  Ll*l » IHH 

IPiT0MAX-TO2PlLl61Illfi9. 190.190 

189  T0MAX«T02P«Litll 

LM-L161 

190  CONTINUE 

PP1NT  2020.T0MAX.LM 
T0MJN»T02PI1) 

LM»i 

1*1  00  199  LI-I.Ihn 

IP<T0MlN-T02P{Ll61 J >199. 199.192 
192  T0MIN«T02P(LUI) 

LM-L161 
199  CONTINUE 

PRINT  2032.TDMIN.LH. 

TWMAX-TWBIl) 

LM-1 

194  00  196  Ll-l.lHM 

IP ( T4MAX-TW8 ( L1L1 ) 1195. 196. 196 

195  TWMAX-TWBI L161 ) 

LM-L161 

196  CONTINUE 

PRINT  2021.TWMAX.LM 
TWMIN*TWB( 1) 

LM-1 

197  00  199  H-I.IMM 

...  IP< TwMIN-TwB IL161 ) 1 199. 199» 198 
196  TWMlN*rwe<L161 I 
LM'ULl 

199  CONTINUE 

PRINT  2039.TWMIN.LM 
ETMAX-TEFFll  J 
LM-1 

200  00  20<  Ll-l.lHM 

IP IETNAX-TCFF <  Lltl I >201 .202 .202 

201  ETMAX-TEFFtLlU  > 

CM-LU1 

202  CONTINUE 

PRINT  2022.ETNAX.LM 

etmin-tcff ( l I 

LM-1 

20J  BO  205  Ll-l.lHH 

1 »  1 1 T M 1  A*—  TtFF"  (Ll^l  1  >205.205*204 

204  ETMIN-TCFF  (LIU  I 
LM«L161 

205  CONTINUE 

PRINT  2094.ETMIN.LM 


tt-U 


V/  \J 


*• 


PRINT  2023 
PRINT  2024 
TDW-0.0 
TWW-0«0 
TEW'0.0 

00  213  J*1 * 1MRS 
T0W«TDirf6TD2P( JI 
TWW>TWW6TWB(JI 
213  T£W»TEW6TEFFO) 

TOS-TOW/MRS 

TW$»TWW/HRS 

TE$«TEW/HRS 

PRINT  2035.TD5, #TwS»TES 
1 P ( NODAYS  J  5000  #  209  #  206 
206  10*1 
L£«24 

00  2071  I0AY*1«N00AYS 
00  207  J-l.fl.LE 
T00(I0AV)*TDDII0AY)6T02PU» 

TWO!  I OAY ) "TmD ( I  DAY )  6TWB ( J ) 

207  TEOIIOAYJ-TEOUOAYI&TEFFIJ) 

LB-LBG24 

2071  LE-LE624 

PRINT  2026 

DO  206  IDAY-l.NOOAYS 
TDAVGI JDAY l*TOO I !DAYl/24.0_ 

TWAVGI I0AY>*TW0( IOAY> /24«0 
TEAVGI  IO/.Y  J  «*TED<  I  DAY  I  /24.0 

208  PRINT  2026 • I OAY *  TOAVGi I OAY ) i TWAVo t I OAY ) i TEAVGI l DAY  > 

209  IFINOINT) 6000 iSOOO *210 

210  PRINT  2027 
NO»2«NOINT“1 

00  212  NI • 1 • ND « 2 
MH* (Nltl 1/2 
ILT-INTINI I 
ILLL* InT I N16 1 ) 

00  211  MI-lLltILLL 
T0I<MMl-TD!(MMl6T02OINI > 

TWI IMMI *TWl I  HMIGTWB IMl I 
'll  TCI I  MM ) *TE1 ( MM ) 6TEFP ( MI I 

T0IAVINM)*T01«MMJ/ I  INTIN  161 l-INTIM  >U > 

TwIAV(HM)>7wI IMKJ/C IN T « N 1 6  *  I - 1  NT ( N I  1 6 1 > 

TCIAV(NMJ«TEI (MMJ/I InTINIGI l-INTIM >61 > 

212  PRINT  202«»INT(NI)  •  INKNIU  I  .TOIAVIMHJ  .T4IAWIMHI  .TEiAVCMMl 
IFUJJJ-J0ATAI4J  1186. 5000*5000 
166  JJJJ* JJ.'Jlr  1 
GO  TO  1311 


1001  FORMAT (F7»JiF;«1 I 

100)  FORMAT (F7*1.F4«1.F5«1.I2.?F7<1*F5.0.F9*2*2F6*1i2F3*0*2I2I 
1004  FORMATUJDI 

1007  FORMAT  IF  5.2 1 

1008  FORMAT  t F4« 1 ) 

1009  FORMAT(2F6«3»2F6.2iF4. JI 

1013  F0RNAT(F9.2.F4,l.F8.3.F4,2.F9*2»F4.l.F4.2»2F4.11 
1016  FORMAT (F5«0*2F9«2*<F6»1I 

1019  FORMAT I2F4«1 *F)»2I 

1020  FORMAT (F7«1i2I3*2F5<2.F)«0> 
lr.21  rORMATIAFA.l  *2F5*2> 


2001  FORMAT  ) lH  ««R  DURING  H0L'R*#13»*  IS  *»F7.4#*  WHICH  is  LESS  THAN  1*0 

l*  therefore#  the  latent  energy  introduced  by  the  occupants  has  bee 

2N  REDUCED**) 

2002  FORMAT) lHl#25X. ‘NOMINAL  OCCUPANT  ACTIVITY  LEVEL*»7X»*400  BTU/HR 
1-OLCUPANT*) 

2003  FORMAT ( lHO *68.v  #*LOADS  (BTu/HR)*) 

2004  FORMAT ( 1kO»7X**V£NT*#59X**TRAN$*  METABOLIC  METABOLIC*) 

2003  FORMAT ) 1H  »*1 IME  RATE  DBT  WBT  CT  AIRCOND  BOUNDARY 
1  CONDEN,  SOLAR  SENSIBLE  LATENT  LIGHTING  EQUIPMENT 
2*1 

2006  FORMAT ( 1H  tt4#3X*4(F4«l#3X)»8(F6*0t3X)) 

2007  FORMAT  I lH  **1NCL I NATION  ANGLE  XI ( •» 12 •• )  IS  «*F6*2#»*  THIS  VALUE  M 
1UST  BE  90  DEG*  OR  LESS**) 

200B  FORMAT ( 1H  **ThE  LATITUDE  OF  THE  SHELTER  IS  GREATER  THAN  THE  PROGRA 
1M  CAN  HANDLE  FOR  THIS  OAY  OF  THE  YEAR**) 

2010  FORMAT  I lH  »*THE  LIMIT  ON  LAYMAX  OR  ON  NMAX  OF  BOUNDARY  **I2»*  HAS 
1BEEN  EXCEEDED**) 

2011  FORM'.TUH  **THE  NUMBER  OF  SHELTER  BOUNDARIES#  NW#  HAS  BEEN  EXCEEDS 
ID#  MACHINE  ERROR**) 

2012  FORMATllH  »*SHELTER  EFFECTIVE  TEMPERATURE  DURING  H0UR*»I4»*  IS  «F5 
1*1**  WHICH  EXCEEDS  THE  LIMITS  OF  THE  METABOLIC  RELATIONS**) 

2014  FORMATllH  »*0RY“8UL8  TEMPERATURE  DURING  HOUR  *»I3«*  IS  **FS*1»*  WH 
1 ICH  EXCEEDS  The  LIMITS  OF  THE  METABOLIC  RELATIONS**) 

2013  FORMATllH  *•  THE  NUMBER  OF  ITERATIONS  FOR  AN  INCREMENT  HAS  EXCEEDS 
ID  The  VALUE  OF  I  DATA! 3) •*) 

2016  FORMATllH  ##TME  RELATIVE  HUMIDITY  IN  THE  SHELTER  HAS  EXCEEDED  SATU 

l RATION#*)  . 

2017  FORMATllH  #*THE  NUMBER  OF  INCREMENTS  IN  THE  HOUR  EXCEEDS  THE  VALUE 
1  OF  I  DATA ( 5 ) * ) 

2015  FORMAT) 1M  *«THE  HOUR  NUMBERS  ARE  NOT  INDEXING  PROPERLY.*) 

2019  FORMATllH  #*THE  HOUR  NUMBER  HAS  EXCEEDED  THE  LENGTH  OF  OCCUPANCY** 
1) 

2020  FORMAT) 1H0#15X#*MAX.  0RY-8ULB  TEMP.  DURING  OCCUPANCY* *6X #F4. 1 .*  AT 
1*#I4.*  HOURS  AFTER  BEGINNING  OF  OCCUPANCY*) 

2021  FORMAT « 1M0#1 3.  ••MAX.  WET-BULB  TEMP.  DURING  OCCUPANCY* *6X *F4. 1 . *  AT 
1 •• 14 1*  HOURS  A*TER  BEGINNING  OF  OCCUPANCY*) 

2022  FORMATUHCISX.OMAX  EFFECTIVE  TEMP.  DURING  OCCUPANCY*  *6X  *F4,  \  .*  AT 
l*#14«*  HOURS  AFTER  BEGINNING  OF  OCCUPANCY*) 

2023  F0RMATHH1.39X.*T1ME-AVERAGE  TEMPERATURE  VALUES*) 

2024  FORMAT»lM0.26X.*FOR  ENTIRE  OURATION  OF  OCCUPANCY*. 16Xi*03T*.7X.*wB 
1T*.»X.*ET* J 

2023  FORMAT ( lHC .26X  **FOR  EACH  DAY  OF  OCCUPANCY*) 

2026  FORMATllH  .6AX.*DAY  * • 12 » 3X »F5 . 1 #2! 5X *F5. 1 ) ) 

2027  FQRMATllMOi2CX.*FOR  SPECIFIC  INTERVALS*»16X»*HRS*) 

2026  FORMATUH  *39X . I  4.*-* . 14 *5X # F3. 1 #2 < 5X #F3. I ) ) 

2030  FORMATUH  .*  HR  CFM/  OEG  OEG  OEG*) 

2031  FORMAT ) lH  «8X  #*OCC  F  F  F*l 

2032  FORMAT) 1M0 • 1 3* **MIN.  ORY-B'/LO  TEMP.  DURING  OCCUPANCY* .6XiF4.li*  AT 
i**I4»*  HOURS  AFTER  BEGINNING  OF  OCCUPANCY*) 

2033  FORMaT(1mO*15X»*MIN.  w£T-  BULB  TEMP.  DURING  OCCUPANCY* .6X «F4 . 1  .*  AT 
1 *. 1 4 1 •  hours  after  beginning  of  OCCUPANCY*) 

2034  FORMAT )1H0«13X»*M1N  EFFECTIVE  T|MP.  DURING  OCCUPANCY*  »6>,  *FA.  1  .*  Af 
1**14.*  HOURS  ArtCR  BEGINNING  OF  OCCUPANCY*! 

2033  FOr.MATtlM  <  7|X  *F3*  1  *2 1 3X.F3. 1  > ) 

2036  FORMAT)lH  **T02P«**F10.1»*  AT  *.l|) 

2037  FORMAT) 1H0.36X.*VALUCS  ARC  FOR  LAST  TIME  INCREMENT  OF  C'CM  HOUR*) 

|03«  FORMAT  UX) 
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